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1  Introduction 


Classical  Gaussian  quadrature  rules  are  extremely  efficient  when  the  functions  to  be 
integrated  are  well  approximated  by  polynomials.  When  the  functions  to  be  integrated 
are  very  different  from  polynomials,  Gaussian  quadratures  do  not  perform  well;  a  par¬ 
ticularly  difficult  problem  involves  the  integration  of  functions  of  the  form 

=  (1) 

t=i 

where  each  of  the  fvmctions  (pi  has  its  own  singularity  at  one  of  the  ends  of  the  interval, 
and  the  function  /  can  only  be  evaluated  tn  toto,  the  coefficients  a,-  being  xmavailable. 
This  problem  is  encountered  in  the  solution  of  integral  equations  with  singular  kernels, 
in  the  numerical  complex  analysis,  in  the  numerical  solution  of  elliptic  partial  differential 
equations  on  regions  with  comers,  and  in  many  other  situations.  While  such  problems 
are  normally  dealt  with  by  means  of  various  ad  hoc  procedures  (see,  for  example,  [1], 
[7]),  these  schemes  lack  the  rapid  convergence,  stability,  and  elegance  of  the  Gaussian 
rules. 

In  fact,  in  [6],  a  far-reaching  generalization  of  the  classical  Gaussian  quadratures  is 
introduced,  replacing  the  polynomials  with  functions  from  an  extremely  wide  class.  The 
quadrature  rules  of  [6]  possess  most  of  the  desirable  properties  of  the  classical  Gaussian 
integration  formulae,  such  as  positivity  of  the  weights,  rapid  convergence,  mathematical 
elegance,  etc.  Unfortimately,  it  is  not  clear  from  [6]  how  such  quadrature  rules  can  be 
obtdned  numerically. 

In  this  paper,  we  present  a  numerical  scheme  for  the  constmction  of  such  generalized 
Gaussian  quadratures.  The  algorithm  is  applicable  to  a  variety  of  functions,  including 
smooth  fvmctions  (not  necessarily  polynomials),  as  well  as  functions  with  end-point 
singularities. 

The  paper  is  organized  as  follows.  In  Section  2,  we  summarize  the  relevant  results 
from  [6],  and  in  Section  3,  we  restate  some  numerical  methods  to  be  used  in  this  paper. 
In  Section  4,  we  develop  analytical  apparatus  to  be  used  in  the  numerical  construction 
of  the  generalized  Gaussian  quadrature  rules,  then  we  extend  those  analytical  tools  to 
functions  with  end-point  singularities  in  Section  5.  The  actual  numerical  algorithm 
is  presented  in  Section  6,  and  the  performance  of  the  algorithm  is  demonstrated  with 
numerical  examples  in  Section  7. 

2  Mathematical  Preliminaries 

In  this  section,  we  summarize  several  classical  results  from  [6]  to  be  used  in  Sections  4 
and  5. 


2.1  Chebychev  Systems 


Definition  2.1  (Ckebyskev  System) 

A  finite  sequence  of  functions  will  be  referred  to  as  a  Chebyshev 

system  if  and  only  if  each  of  them  is  continuous  on  [a,  6],  and  the  determinants 

V’lM  •••  y>i(®m) 

(Piixi)  (fia{x2)  •••  (pilxn) 

*  »  • 

•  •  • 

•  •  • 

'Pmixi)  (fmiXi)  •••  (Pm{Xm) 

are  non-zero  for  any  set  of  m  points  x\,X2,  -  •  •  tXm  €  [o,  6]  such  that  x,-  ^  Xj  for  any 
ii^j. 

Following  are  several  important  cases  of  Chebyshev  systems  (for  more  examples,  see 

[6]). 

Example  2.1  For  any  natural  m,  the  functions  1,  x,  x^,  •  •  • ,  x”*  constitute  a  Chebyshev 
system.  Moreover,  if  oi\,  012,  ol^  is  a  sequence  of  distinct  real  numbers,  then  the 
system  {x^**}  is  a  Chebyshev  system  on  any  interval  [a,  6]  C  (0,oo). 

Example  2.2  For  any  n  distinct  real  numbers  oriyas)  •  *  *  the  functions  c"**,  xe“‘*, 
c®’*,  xc®’*,  •  •  -,  c®"*,  xc®"*  constitute  a  Chebyshev  system  on  any  interval  [o,  6]  C 
(—00,00). 

Definition  2.2  (Hermite  System) 

A  finite  sequence  of  functions  be  referred  to  as  an  Hermite 

system  on  the  interval  [a,  b]  if  and  only  if  (pi  €  C^fa,  6]  /or  a// 1  =  1, 2,  ■  •  • ,  2n,  and  the 
determinants 

V\{X2)  <Pi(x2)  •••  V>x(®n)  V»i(afn)  ^ 

<P2ixi)  /j(xi)  (P2{X2)  q>2{.^2)  •  •  •  q>2{xn)  V2M 

:  :  :  :  :  : 

¥»2n(ari)  V’anCaJl)  V’2n(a;2)  ^2n{^2)  ■”  V>2n(Xn)  V»3„(x„)  > 

ore  non-zero  for  any  set  of2n  points  xi,X2,- •  •  ,X2n  €  [a, 6]  such  that  x,  ^  xj  for  any 

i¥^j- 

Definition  2.3  (Extended  Hermite  System) 

A  finite  sequence  of  functions  {¥>i,V’2j"-»V’2n}  will  be  referred  to  as  an  extended 
Hermite  system  if  it  is  both  Chebyshev  and  Hermite. 

Remark  2.1  The  Extended  Hermite  systems  are  a  slight  generalization  of  the  extended 
Chebychev  systems  of  [6]. 

Following  are  several  important  cases  of  extended  Hermite  systems  (for  more  exam¬ 
ples,  see  [6]). 
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Example  2.3  For  any  natural  n,  the  functions  constitute  an  extended 

Hermite  system.  Moreover,  if  a\,  a2,  •••,  oc^n  is  a  sequence  of  distinct  real  numbers, 
then  the  system 

•  •  ,x“*"  (4) 

is  an  extended  Hermite  system  on  any  interval  [a,  6]  C  (0,  cx)). 


Example  2.4  An  important  special  case  of  the  proceeding  example  is  the  finite  sequence 
of  functions 


l,x^x,x^+^x^x"+^■ 


n-1  n-l+or 

>  *  j  * 


(5) 


with  a  an  arbitrary  non-integer  real  number. 


2.2  Gaussian  Quadratures 

We  will  be  considering  integrals  of  the  form 

I  oj(x)<p{x)dx,  (6) 

Ja 

where  u)  :  [a,  b]  —*  R}  is  a.  non-negative  function  to  be  referred  to  as  the  weight  function, 
and  (fi  :  [a,  6]  R}  is  &  function  from  a  suitably  chosen  class.  A  quadrature  rule  is  an 
expression  of  the  form 

Tn{^)  =  ^Wi-<p{Xi),  (7) 

»=1 

with  X,-  €  [a,  6]  and  w,  €  R}  for  all  t  =  1, 2,  *  •  * ,  n.  The  points  x,  and  coefficients  Wi  are 
referred  to  as  the  nodes  and  waghts  of  the  quadrature  formula  (7),  respectively,  while 
the  expression  (7)  itself  is  viewed  as  an  approximation  to  the  integral  (6).  Normally, 
quadrature  formulae  are  chosen  to  be  exact  on  certain  chosen  sets  of  functions,  most 
frequently,  polynomials  up  to  some  fixed  order  m.  An  n— point  quadrature  formula  is 
referred  to  as  a  Gaussian  quadrature  if  and  only  if  it  integrates  exactly  all  polynomials 
of  orders  up  to  2n  —  1. 

We  will  generalize  the  notion  of  the  classical  Gaussian  quadrature  somewhat,  by 
introducing  the  following  definition. 

Definition  2.4  (^Gaussian  Quadrature) 

Suppose  that 

• '  ■>92n}  (8) 

is  a  set  of  integrable  functions  [a,  h]—*R}.  We  will  say  that  the  n— point  quadrature  rule 
(7)  is  Gaussian  with  respect  to  the  system  (8)  if  and  only  if  it  integrates  exactly  all  of  the 
functions  (8).  In  other  words,  a  Gaussian  rule  is  an  n— point  rule  that  is  exact  for  2n 
pre-chosen  functions.  We  will  refer  to  the  nodes  and  weights  of  a  Gaussian  quadrature 
as  the  Gaussian  nodes  and  weights,  respectively. 
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Remark  2.2  Obviously,  a  classtad  Gaussian  quadrature  rule  is  a  Gaussian  quadrature 
rule  for  which 

<pi(x)  =  1, 

IP2(X)  =  X, 


(fi2n{x)  =  X^“‘. 


(9) 


The  principal  result  we  use  from  [6]  is  the  following  theorem. 


Theorem  2.1  (KarUn-Studden) 

Suppose  that  the  functions  92,  *  ■  ■ ,  constitute  a  Chebyshev  system  on  the  inter¬ 
val  [a,  6].  Then  there  exists  a  unique  n-point  quadrature  rule  (7)  that  is  Gaussian  with 
respect  to  the  functions  {91,  ^3,  ■  ** ,  ^3n}-  Furthermore,  all  the  weights  wi,  103,  ■  ■  ■ ,  u’n  of 
the  quadrature  are  positive. 


As  for  smooth  functions,  Gaussian  quadrature  rules  also  exist  for  a  variety  of  func¬ 
tions  with  end-point  singularities.  The  following  theorem  is  an  immediate  consequence 
of  Theorem  2.1. 


Theorem  2.2  Suppose  that  functions  :  (a,  6]  — »  Ff  are  continuous,  and  integrahle 
on  [a,  6]  for  a// 1  =  1,2,  ■  *  •  ,2n.  Suppose  also  that  the  function  w(x)  >  0  is  continuous 
on  (a,  6]  and  integrahle  on  [a,  6].  Suppose  further  that  functions  t/>i  are  defined  by  the 
formula 


i/>i{x) 


(10) 


and  that 


lim  rfftix)  <  00 

►a  '  ' 


(11) 


for  all  i  =  1,2,  •••,2n.  Suppose  finally  that  the  functions  »02n}  defined  by 

(10)  constitute  a  Chebyshev  system  on  the  closed  interval  [a,  6]. 

Then  there  exists  a  unique  n-point  quadrature  rule  (7)  that  is  Gaussian  with  respect 
to  the  functions  {<pi,<p3,"  ’  ,<P2n}-  Furthermore,  all  the  weights  wi,W2,' "  ,von  of  the 
quadrature  are  positive. 


Proof:  The  theorem  is  proved  by  applying  Theorem  2.1  to  the  new  weight  function 


<i>(x)  =  ii){x)  •  u;(z), 

and  the  new  set  of  functions  {^1,  ^2i'"i  ■ 


(12) 


Example  2.5  For  any  natural  n,  and  real  number  0  <  a  <  1,  the  unique  n-point 
quadrature  on  the  interval  [0, 1]  with  respect  to  the  functions 

1,  X,  x*“*,  X*,  x*“®,  •  •  • ,  x**"^ ,  x""'"®. 

can  be  obtained  via  the  following  Chebyshev  system  (see  Example  2.3) 

X®,  1,  X*+®,  X,  x*+®,  x^ . . . ,  x"-^+®,  x"-^ 

on  the  interval  [0, 1]  with  the  weight  function  u;(x)  =  u(x)  •  x®. 


(13) 

(14) 
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3  Numerical  Preliminaries 


In  this  section,  we  collect  the  relevant  ntunerical  tools  to  be  used  in  Sections  4  and  5. 
They  can  be  found,  for  example,  in  [3],  [4],  [5]. 

3.1  Nested  Chebyshev  Approximation 

For  any  non-negative  integer  n,  the  Chebyshev  polynomial  Tn  of  order  n  is  defined  by 
the  formula 

Tn{cos0)  =  cos{n0).  (15) 

Clearly,  |rn{x)|  <  1  for  z  €  f-1,1]. 

The  Chebyshev  polynomials  constitute  an  orthonormal  basis  for  I>^[— 1,1]  with  re¬ 
spect  to  the  inner  product 

{f,9)  =  ‘ 9{x)dx.  (16) 

Therefore,  any  function  /  €  C®[— 1, 1]  can  be  represented  by  an  expansion 

/(x)  =  f:<n-Ti(i),  (17) 

»=0 

with  the  coefficients  Oi  pven  by 

<U  =  (/,?()•  (18) 

Lemma  3.1  states  that  the  Chebyshev  series  (17)  converges  rapidly  for  sufficiently 
smooth  functions.  Its  proof  can  be  found,  for  example,  in  [5]. 

Lemma  3.1  Suppose  that  n  and  k  are  natural  numbers,  and  that  f  €  1).  Sup¬ 

pose  further  that  the  coefficients  ao,oi,--*,a„  are  defined  the  formula  (18).  Then  for 
any  X  €  [-1,1], 

m-±atT,{x)  (19) 

vsO  ” 

In  particular,  if  f  &  C°°,  then  the  expansion  (17)  converges  to  f  superalgebraically. 

Observation  3.1  For  functions  uhth  end-point  singularities,  such  as  /(z)  =  Inz,  we 
can  build  a  structure  on  the  given  interval,  consisting  of  subintervals  clustering  near  the 
end  points  (see  Figure  3.1),  and  then  use  the  Chebyshev  expansion  (17)  to  approximate 
the  functions  on  each  subinterval.  On  each  of  the  subintervals,  the  Chebyshev  expansion 
converges  superalgebraiclly. 
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Figure  1:  Subiutervals  clustering  near  the  end-points 

3.2  Nonlinear  Equations 

We  will  be  considering  systems  of  nonlinear  equations  of  the  form 

F{x)  =  0, 

where  x  =  (®i,  •  •  • , Xn)^  €  i2“,  and  the  mapping  F  :  iP*  — »  fT*  b  of  the  form 


F(x)  = 


flixir'-iXn) 
hiXU‘"  tXn) 


V  “  *  *  >  ®n)  / 

Definition  3.1  The  Jacobian  matrix  of  mapping  F  in  (21)  is  defined  by  the  formula 

/  fh.  ...  |ii  \ 

S*1  Sxn 

DF(x)  =  ;  :  .  (22) 

,  iJa.  . . .  fik 

'  9zi  dSn  ' 

The  following  wo  lemmas  about  the  solution  of  systems  of  nonlinear  equations  (20) 
are  well-known  ^see  [3],  [8],  for  example). 

Lemma  3.2  (Newton’s  Method) 

Suppose  that  F  :  IV*  IP*  is  continuously  differentiable  in  an  open  convex  set 
D  C  IC*,  and  the  mapping  G  :  IV*  -*  IP*  is  defined  by  the  formula 

y  =  X  -  (DF(x))-'F(x).  (23) 

Suppose  also  that  x*  £  IV*  is  the  zero  of  F,  and  there  exists  >  0  such  that 

||i?F(xrMl  < /».  (24) 

Suppose  further  that  there  exist  two  positive  numbers  r  and  7  such  that  x  €  Z)  for  any  x 
such  that  ||x  —  x*||  <  r,  and 

||DF(x)-DF(y)||<7l|x-yl|  (25) 

for  any  x,  y  such  that  (|x— x*||  <  r,  |(y  —  x*((  <  r.  Suppose  finally  that  xq  is  an  arbitrary 
point  in  IV*,  and  the  sequence  Xi,X2,  •  •  of  points  in  IV*  is  defined  by  the  formula 


Xfc+i  =  G{xk) 
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for  allk  —  0, 1,2, ••  •• 

Then  there  exist  e  >  0  and  a  >  0  such  that  the  sequence  generated  by  (26)  converges 
to  X*,  and 

llxk+i  -  x’ll  <  a||xk  -  x*||*  (27) 

for  any  Xq  such  that  ||xo  —  x*||  <  c. 

Lemma  3.3  (Modified  Newton’s  Method) 

Suppose  that  under  the  assumptions  of  Lemma  S.2,  x*  is  the  zero  of  the  mapping 
F  :  D  —*  Suppose  also  that 

AQ,Ai,A2t-“  (28) 

is  a  sequence  ofnxn  nonsingular  matrices,  and  the  mapping  G  :  —*  R*  is  defined 

by  the  formula 

y  =  x-A;^F(x).  (29) 

Suppose  further  that  there  exists  a  positive  real  number  M  such  that 

iiA.-Df(xi)||<M-||P(xk)||.  (30) 

Suppose  finally  that  Xq  is  an  arbitrary  point  in  and  the  sequence  Xi,X2,  •  *  •,  of  points 
in  RP'  is  defined  by  the  formula 

x*+i  =  G{xk)  (31) 

for  a//  A:  =  0, 1,2,  •  • 

Then  there  exist  e>  0  and  a  >  0  such  that  the  sequence  generated  by  (SI)  converges 
to  X*,  and 

I  |xk+i  -  X*  1 1  <  a|  |xk  ~  X*  1 1*  (32) 

for  any  Xo  such  that  ||xb  —  x*||  <  t, 

3.3  Continuation  Method 

The  Newton  algorithm  for  the  solution  of  systems  of  non-linear  equations  is  an  extremely 
powerful  technique,  provided  that  a  satisfactory  initial  point  is  available.  In  many  cases, 
a  starting  point  is  not  available  directly,  but  can  be  obtained  by  the  process  known  as 
the  continuation  method  (otherwise  referred  to  as  the  homotopy  method).  Following  is 
a  brief  description  of  the  technique. 

Suppose  that  we  are  trying  to  solve  a  system  of  non-linear  equations 

F(x)  =  0,  (33) 

with  F  :  BP  BP  satisfying  the  conditions  of  Lemma  3.3,  except  for  the  initial  point 
Xo,  which  is  not  available.  Suppose  further  that  we  do  have  2u:cess  to  a  mapping 
G  :  [0, 1]  X  ^  ,  satisfying  the  following  conditions. 
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1.  The  mapping  Fo  :  if*  — »  if'  defined  by  the  formula 

Fo(x)  =  G(0,x)  (34) 

has  a  simple  structure,  so  that  the  solution  of  the  equation 

Fo(x)  =  0  (35) 

is  unique  and  known. 

2.  For  all  X  €  R^, 

G(l,x)  =  F(x).  (36) 

3.  For  all  t  €  [0, 1],  the  equation 

C?(t,x)  =  0  (37) 

has  a  unique  solution  x(t)  €  F",  and  satisfies  the  conditions  of  Lemma  3.3  in  the 
neighborhood  of  x(t). 

4.  X  is  a  continuous  (or  better,  Lipschitz)  function  of  t. 

Under  the  above  conditions,  the  following  procedure  yields  the  solution  of  the  equa¬ 
tion  (33). 

1.  For  a  sufficiently  large  m,  construct  the  points  ti  =  (t  —  l)/(m  —  1)  on  the  interval 
(0, 1],  with  i  =  1, 2,  •  •  • ,  m,  and  consider  the  solutions  of  the  equation  (33)  for  <  =  t < 
with  j  =  1,2, •  •  •  ,m. 

2.  Clearly,  we  know  the  solution  x(0)  of  the  equation, 

(?(ti,x)  =  0  (38) 

when  *  =  1,  and  for  all  *  =  2,3, •  •  •  ,m,  we  solve  the  equation  (37)  by  means  of 
Lemma  3.3,  using  x(t,--i)  as  the  initial  approximation. 

3.  Since  =  1,  the  result  of  the  final  step  of  this  process  is  the  solution  of  the 
equation  (33). 

A  detailed  discussion  of  the  continuation  techniques  can  be  found  in  [4],  where  the 
convergence  of  the  above  scheme  is  proven  (in  a  much  more  general  environment)  for  all 
stifficiently  large  m. 

4  Analytical  Apparatus 

In  this  section,  we  develop  analytical  tools  to  be  used  in  the  numerical  construction  of 
the  Gaussian  quadratures  whose  existence  follows  from  Theorem  2.1. 
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4.1  Construction  of  Gaussian  Quadratures 

The  following  lemma  is  an  immediate  consequence  of  Definition  2.3  of  extended  Hermite 
systems. 

Lemma  4.1  Suppose  that  the  functions  {(fi,  (P2t  *  ■  * » fin}  constitute  an  Hermite  system 
on  the  interval  [n, 6]^  und  xi,X2, •  •  •  ,a;n  ^  points  on  the  interval  [a, 6]  such  that 
Xi  ^  Xj  for  any  i^  j.  Then  there  exist  such  unique  coefficients  Qij,  Aj,  t  =  1, 2,  •  •  • ,  n, 
j  =  1,2, •  •  •  ,2n  that 

f  ».;(!»)  =  0,  (39) 

<T,(X*)  =  Oik, 

I  ^  n*’  (40) 

I  ViM  =  0, 

/or  a// 1  =  1, 2,  •  • ,  n,  ib  =  1, 2,  •  •  • ,  n,  with  Si,k  denoting  the  Kronecker  symbol,  and  the 
functions  ai,  rfi  defined  by  the  formulae 

’  ^■'(®)’ 


Viix)  -  £  Aj  •  ‘Pjix)-  (42) 

ial 

Furthermore,  there  exist  unique  coefficients  a^,  bij  with  i  =  1, 2,  •  •  * ,  2n,  j  =  1, 2,  •  •  • ,  n, 
such  that  ^ 

j=i 

for  all  i  —  1, 2,  •  •  • ,  2n. 

In  other  words,  there  exist  unique  linear  combinations  o’i,o’2, •  •  •  ,cr„,  qi,V2,'  ’  •  ,Vn 
of  the  functions  <Pi,<P2,’ "  ,V2n  satisfying  the  conditions  (39),  (40),  (40)-  Conversely, 
the  functions  <^1,(^2, '**  7*^2n  ore  linear  combinations  of  the  functions  (ri,o’2, •  •  •  ,0’n> 

Vi,V2,--'  ,Vn- 

Theorem  4.1  below  is  the  principal  analytical  tool  of  this  paper.  It  establishes  the 
necessary  and  sufficient  conditions  und^  which  a  quadrature  is  Gaussian  with  respect 
to  a  given  Hermite  system. 

Theorem  4.1  Suppose  that  functions 

<P\,<P2,'”'>V>2n  (44) 

constitute  an  Hermite  system  on  the  interval  [o, &].  Then  the  nodes  Xi,X2, •  •  ■  jXn  on 
[a,  6]  are  Gaussian  with  respect  to  the  functions  (44)  */  ond  only  if 


fh 

J  (jj{x)ai{x)dx  =  0 
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for  alii  =  1,2,’  ’•  ,n.  In  this  case,  the  Gaussian  weights  xo\,Wi,  *  *  * ,  tOn  are  given  by  the 
formula 

Wi  =  J  uj{x)rii{x)dx,  (46) 

with  functions  (Tj,  •  •  • ,  <t„,  »yi,  172,  •  •  *  >  Vn,  *«  (45),  (45)  are  defined  by  formulae  (39)  - 
(42). 


Proof:  First,  we  show  that  for  any  Gaussian  quadrature  rule  (7),  the  conditions  (45), 
(46)  are  satisfied.  Indeed,  since  it  integrates  exactly  all  of  the  functions  ••• ,  (^2n) 

it  also  integrates  exactly  all  their  linear  combinations  <ri, <r2,  “  -  >  i7i>  ’?2)  *  *  *  >  Vn-  Now, 
(45),  (46)  follow  immediately  from  (39),  (40). 

Suppose  now  that  the  nodes  xi,X2, are  such  that  the  conditions  (45)  are 
satisfied,  and  the  coefficients  u}i,Wi,.  •  •  ,Wn  are  defined  by  the  formula  (46).  We  will 
show  that  the  n-point  quadrature  (7)  is  Gaussian  with  respect  to  the  system  (44). 

Due  to  Lemma  4.1,  there  exist  coefficients  Qij,  0ij  ,  i  =  1, 2,  •  •  • ,  n,  j  =  1, 2,  •  •  • ,  2n 
such  that 

n 

i=i 

for  any  t  =  1, 2,  •  •  • ,  2n.  Thus 

fb  rb  rb 

j  u{x)<pi{x)dx  —  ^  a,j  /  u;(x)<T,(x)dx  +  6,,  /  u{x)r)j{x)dx.  (48) 
Ja  Ja  Ja 

Combing  (48)  with  (45),  (46),  we  have 

/  u>(x)v?,(x)dx  =  22  (49) 

i=i 

On  the  other  hand,  combining  (39),  (40),  and  (48),  we  obtain 

n  n  n 

S  H  “'i  2(a.*o-*(x,)  +  bikVkixj) 

3=1  fc=i 

=  (50) 

i=i 


Combining  (49)  and  (50),  we  finally  get 


/ 


u}{x)(fi{x)dx  = 


(51) 


for  all  t  =  1,2, •  •  •  ,2n.  ■ 

Theorem  4.2  below  follows  immediately  from  Theorem  4.1.  It  describes  the  Gaussian 
nodes  as  the  solution  of  a  system  of  non-linear  equations. 
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Theorem  4.2  Suppose  that  the  functions  •  •  •  ,yjn}  constitute  an  Hermite  sys¬ 

tem,  and  functions  {<ri,  <T2,  •  •  • ,  <j„}  and  {rii,T]2,  •  •  • ,  i/n}  defined  by  the  formulae  (89), 
(40),  (4l)>  (48)-  Suppose  further  that  S  is  a  subset  of  R*  consisting  of  all  finite 

sequences  {xijXjj '  •  *  such  that  Xi  ^  Xj  whenever  i  ^  j  .  Suppose  finally  that  the 
mapping  F  :  S  is  defined  by  the  formula 


^  Xi  ^ 

^  Xfu;(x)o'i(x)dx  ^ 

X2 

— 

fj^u{x)or2{x)dx 

V  a:„  J 

<  Sa  (^{x)(Tnix)dx  j 

(52) 


Then  {xi,X2,  •  •  •  ,2n}  ere  the  Gaussian  nodes  with  respect  to  the  system  of  functions 
{(pi,<p2,  ■  •  • ,  y>n}  if  and  only  if 

F{xi,---,Xn)  =  Q.  (53) 


4.2  Computation  of  Gaussian  Quadratures 

In  t'.'s  subsection,  we  observe  that  the  modified  Newton  method  in  Subsection  3.2 
assumes  a  particularly  simple  form  when  it  is  applied  to  the  system  of  equations  (53), 
and  establish  Theorem  4.5,  the  principal  numerical  tool  of  this  paper.  Theorem  4.5 
shows  that  an  extremely  simple  iterative  scheme  converges  quadraticaUy  for  the  system 
of  equations  (53). 

Theorem  4.4  below  provides  an  analytical  expression  for  the  Jacobian  matrix  of  the 
mapping  F  defined  by  the  formula  (52).  Theorem  4.3  is  the  consequence  of  Lemmas  4.2 
and  4.3,  and  will  be  used  in  the  proof  of  Theorem  4.4.  Theorem  4.5  follows  immediately 
form  Lemma  3.3,  Corollary  4.2,  and  Theorem  4.4. 

Lemma  4.2  Suppose  that  functions  {(pi,(p2,  ,  ¥’2n}  form  an  extended  Hermite  system 

with  (pi  €  C^\a,  6]  for  aW  i  =  1, 2,  •  ■  • ,  2n.  Suppose  also  that 


Xl,X2,"-,Xn  (54) 

are  n  distinct  points  on  the  interval  [a,  6],  and  functions 

0’1>0'2,  (55) 

are  determined  by  the  set  of  points  (54)  via  formulae  (39)  and  (40).  Suppose  further 
that  I  is  an  integer  such  that  I  <  I  <n,  and  6  is  a  real  number  such  that 


+  5,Xf+i  •••,Xn  (56) 

are  n  distinct  points  on  interval  [a,  5].  Suppose  finally  that  the  functions 

5'i,d2,  •••,&„,  7i,  ^2,- (57) 
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art  determintd  by  the  set  of  points  (56)  viaformtdae  (89)  and  (40). 
Then  there  exist  coefficients  a,/  and  with  i  =  1, 2,  *  •  • ,  n,  such  that 


cr<(x)  =  <Ti{x)  +  an  •  <ti{x)  +  fin  •  i7/(x) 

(58) 

for  alli^l,  i  —  1, 2,  •  •  • ,  n,  and 

flr/(x)  =  an  •  <T/(x)  +  fin  •  7/(x). 

(59) 

Proof:  Due  to  Lemma  4.1,  there  exist  unique  coefficients  aij,fiij  with  i 
j  =  1, 2,  ■  •  • ,  n,  such  that  for  all  t  =  1, 2,  •  •  • ,  2n, 

—  1.2, 

n 

•  Vj{x)) 

(60) 

for  any  x  €  [o,5].  Differentiating  (60),  we  have 

^i(®)  =  S  (^i  *  •  O'jix))  . 

(61) 

For  the  functions  (57),  the  conditions  (39),  (40)  assume  the  form 

r  5i{xi  +  6)  =  0, 

\  a'i(xi  +  8)  =  6n, 

(62) 

and 

1  Oi{xk)  =  0, 

\  ollxk)  =  6ik 

(63) 

for  all  t  =  1, 2,  •  •  • ,  n,  fc  =  1, 2,  •  •  • ,  n,  and  k^l. 

For  any  k  ^  1,  evaluating  (60)  at  Xi^  and  applying  conditions  (39), 
obtain 

fiik  =  0 

(40),  (63),  we 

(64) 

for  all  i  =  1, 2,  •  •  • ,  n. 

Similarly,  for  any  k  ^  1,  evaluating  (61)  at  Xk  and  applying  conditions  (39),  (40), 

(63),  we  have 


Oik  =  Sik  (65) 

for  all  *  =  1, 2,  •  •  • ,  n. 

Now,  (58),  (59)  immediately  follow  from  (60),  (64),  (65).  ■ 

Lemma  4.3  Suppose  that  under  the  assumptions  of  Lemma  4.9,  the  coefficients  oiu  and 
fin  are  defined  via  formulae  (58)  and  (59)  for  alii  —  1,2,  •  •  •  ,n,  Then 

0.1  =  + 

A  =  0(«') 
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(66) 

(67) 


(68) 

(69) 


for 


1,2, and 

an  =  l-6-a;{xi)^-0{6^)), 

A/  =  6-\-0{6\ 


Proof:  Expanding  the  functions 

I  i 

^it  Vii  Vii 

into  the  Taylor  series  at  z/,  we  have 


(  (Ti{xi  +  S)  =  <Ti{xi) -h  6  •  <T’i(xi)  +  0{^), 
\  r)i{x,  +  6)  =  ^  +  0(5*), 

/  o<(xi  +  5)  =  a({xi)  +  6 .  <rUx,)  +  0(5*), 
\  i7-(x/  +  5)  =  i;<(x/)  +  5  •  »^"(x/)  +  0(5*), 


for  all  t  =  1, 2,  •  •  • ,  n. 

Evaluating  (58)  and  (59)  at  x  =  (xi  +  5)  and  using  the  conditions  (39),  (40),  (62), 
(70),  we  obtain 


A,  = -a./.  5 +  0(5*)  (72) 

for  all  t  =  1,2, •  •  •  ,n. 

Differentiating  (58)  at  x  =  (xj  +  5)  and  using  conditions  (39),  (40),  (62),  (71),  we 


have 

a„  =  (l+<.<rr(x,))‘'(l  +  0(<’)), 

=  l-S-<T“(x,)-¥0{e‘).  (73) 

Similarly,  differentiating  (59)  at  x  =  (xj  +  5)  and  applying  the  conditions  (39),  (40), 
(62),  and  (71),  we  get 

a, 7  =  -S  •  ct”{x,)  +  0(5*)  (74) 

for  all  i  =  1, 2,  d  •  •  • ,  n,  and  i  ^  / 

Finally,  combining  (72)  with  (73),  (74),  we  have 

fin  =  5  +  0(5*),  (75) 

fin  =  0(5*)  (76) 

for  all  t  =  1, 2,  ’  •  • ,  n,  and  I  M 

Combining  Lemmas  4.2,  4.3,  we  now  obtain  the  following  theorem. 

Theorem  4.3  Under  the  assumptions  of  Lemma  4-2, 

at{x)  =  <T,(x)  -  5  •  (r"(x,)  •  a,(x)  -  5  •  T7j(^)  +  0(5*)  (77) 

=  a,(x)-5-<7"(x,)-(T,(x)  +  0(5^)  (78) 

for  any  t  =  1, 2,  •  •  • ,  n,  and  i^l. 
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The  following  theorem  is  an  immediate  consequence  of  Theorem  4.3.  It  provides  a 
simple  expression  for  the  Jacobian  of  the  mapping  (52),  showing  that  the  latter  is  nearly 
diagonal  in  the  vicinity  of  the  solution  of  the  equation  (53). 

Theorem  4.4  Suppose  that  functions  form  an  extended  Hermite  sys¬ 

tem,  and  <pi  €  C®[a, b]  for  i  =  1,2,  •  •  •  ,2n.  Suppose  further  that  x\,  x^,  •••,  x^  are  n 
distinct  point  on  the  interval  [a,  6],  and  functions 


<3^3>  •  •  *  J  1?3,  *  •  *  ,  »/n 

are  determined  by  formulae  (S9)  and  (40) •  Then  the  Jacobian  DF(x.)  of  the  mapping  F 
defined  by  (52)  is  given  by  the  formula 


DF{x)  =  - 


I  u(x)iji(x)dx  0 

0  J?w(a:)j73(x)<ix 


0 

0 


\  0  0  •••  f^  U}(x)tin(x)dx 

udiere  x  =  (xi,  xj,  •  •  • ,  x„)^,  and  matrix  E(x)  is  given  by  the  formula 


-E(x)  (79) 


/  <ri'(xi)  •••  (Tiixn)\  (  Xfw(x)<ri(x)dx  0 


£;(x)  = 


0 


W«(®i)  ••• 


0 


0  •••  jJa;(x)<Tn(x)dx 


.  (80) 


Froof:  Suppose  that  /  is  an  integer  such  that  I  <  I  <  n,  and  S  is  a  real  number  such 
that 


*!>••*>  +  ^5  Xl+l 


>*n 


(81) 


are  n  distinct  points  on  the  interval  [a,  6].  Suppose  further  that  the  functi(ms 

5’l,5’3,-*-j^n»^l,^3>”'»7»»  (82) 

are  determined  by  the  set  of  points  (81)  via  formulae  (39)  and  (40), 

Combining  (77),  (78)  with  Definition  3.1  of  Jacobian  matrix,  and  the  definition  (52) 
of  the  mapping  F,  we  immediately  obtain 

(DF(x))^i  =  ^  ^  u;(x)<r,(x)dx 


_  /,  u;(x)dt(x)dx  -  /,  u;{x)(Ti{x)dx 

$•^0  S 

=  — /  (jj{x)ai{x)dx  —  /  uj{x)rii{x)dx 
Ja  Ja 


(83) 
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for  any  t  ^  t  =  1, 2,  •  ■  • ,  n,  and 


iDF{x)),, 


~  jf  ij(z)<Ti(x)dx 

Uiu  -  I^(^{x)<ri{x)dx 

8—0  S 

-ffiixt)  j  Ul{x)<Tl{x)dx. 


(84) 


Corollary  4.1  follows  immediately  from  Theorem  2.1  and  Theorem  4.4,  and  corollary 
4.2  is  the  consequence  of  Corollary  4.1  and  Theorem  4.4. 


Corollary  4.1  Suppose  that  under  the  assumptions  of  Theorem  4>4t  function  F  : 

is  defined  by  (52).  Then  there  exists  a  unique  x*  =  (*;,  xj,  •  •  • ,  x*  €  /f* 

such  that 

Fix-)  =  0,  (85) 

and  the  Jacobian  matrix 


DFix-)  =  - 


/  !lf^ix)mix)dx 


0 

!t(^ix)rhix)dx 


0 

Q 


\ 


(86) 


\ 


0 


0  •••  Xfa;(x)f7n(x)dx  ) 


is  nonsingular,  where  the  functions  **  »*?*»  determined  by  the  set  of  points 

xj,  xj,  •  •  • ,  X*  via  the  formula  (40). 


Corollary  4.2  Suppose  that  under  the  assumptions  of  Theorem  4‘4>  function  F  : 

is  defined  by  (52),  and  x*  is  Rie  unique  zero  of  F.  Then  F  is  continuously 
differentiable,  and  there  exist  three  positive  real  numbers  r,  /?  and  7  sudt  that 

||X)F(xV||</S,  (87) 


and 

||Z?F(x)-/>F(y)||<7l|x-y||  (88) 

for  any  x  and  y  such  that  ||x  —  x*||  <  r,  ||y  —  x*||  <  r. 

The  following  theorem  is  the  principal  numerical  tool  of  this  paper.  It  shows  that  an 
extremely  simple  iterative  scheme  converges  quadratically  for  the  S3rstem  of  equations 
(53),  and  is  an  immediate  consequence  of  Lemma  3.3,  Corollary  4.2,  and  Theorem  4.4. 
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Tkeorem  4.5  Suppose  that  functions  form  an  extended  Hermite  sys¬ 
tem,  and  the  mapping  G  :  R*  is  defined  by  the  formula 


with  t  s=  1,2, 


_  f^uf(x)ai(z)dx 


,  n,  and  the  functions 


(89) 


•  1  *?1»  *?2»  •  •  • » ^?n 

defined  by  the  points  xi,  zj,  *  *  ■ ,  Xn  via  the  formulae  (S9)  and  (40).  Suppose  further  that 
<pi  €  C®[a,  b]  for  alii  =  1, 2,  •  •  • ,  2n,  and  the  function  F  :  -*  IC  is  defined  by  (52). 

Suppose  finally  that  is  the  unique  zero  of  F,  that  Xo  is  an  arbitrary  point  in  R*,  and 
the  sequence  Xi,X3,  •  •  •,  of  points  in  is  defined  by  the  formula 

x/h.i  =  G'(x*)  (90) 

for  all  k  =  0,1,2,“  - . 

Then  there  exists  e  >  0  arid  a  >  0  such  that  the  sequence  Xi,X3,  generated  by 
(90)  converges  to  x*,  and 

l|xk+i-x*||<a|jxk-x*||’  (91) 

for  any  initial  point  xo  such  that  ||xo  —  x*{|  <  c. 

Remark  4.1  In  Theorem  4.5,  we  impose  the  condition 

V>,  €  (92) 

for  all  i  =:  l,2,.'‘,2n.  However,  it  can  be  easily  observed  that  the  condition  (92)  is 
excessively  restrictive,  and  a  somewhat  more  involved  proof  shows  that  as  long  as  are 
continuously  differentiable,  and  satisfy  the  Lipschitz  condition  for  a// 1  =  1, 2,  ••  • ,  2n, 
the  modified  Newton  method  (90)  will  stUl  converge  quadratically. 


5  Integration  of  Singular  Functions 

In  this  section,  the  theory  of  the  generalized  Gaussian  quadrature  rules  established  in 
Section  4  wiU  be  generalizeo  a  variety  of  functions  with  end-point  singularities.  We 
will  first  introduce  the  concepts  of  Chebyshev  and  extended  Hermite  systems  in  the 
case  of  singular  functions.  Then  we  wiP  prove  Theorems  5.1  and  Theorem  5.2  for  the 
construction  of  Gaussian  quadrature  ruit;s,  providing  effective  numerical  construction 
for  Gaussian  quadratures  for  functions  with  end-point  singularities. 

Definition  5.1  (Chebyshev  System) 

Suppose  that  X  C  is  either  {a,b),  or  [a, 6),  or  (a, 6].  Then  a  finite  sequence  of 
functions  ¥’2^  *  * '  >  ^n}  vnll  be  referred  to  as  a  Chebysh'iv  system  on  X  if  and  only  if 
it  constitutes  a  Chebyshev  system  on  every  closed  subirAerval  [c,d[C.X  (set  Definition 
2.1). 
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Definition  5.2  (Extended  Hermite  System) 

Suppose  that  X  C  R}  is  either  (a,b),  or  [a,  6),  or  (a,fr].  Then  a  finite  sequence  of 
functions  {^1,92,-  *  ■  ,9n}  toi//  he  referred  to  as  an  extended  Hermite  system  on  X  if  and 
only  if  it  constitutes  an  extended  Hermite  system  on  every  closed  subinterval  [c,d^C  X 
(see  Definition  8.S). 

The  following  is  an  important  example  of  extended  Hermite  systons. 

Example  5.1  For  any  natural  n,  the  functions 

1,  :e“,  X,  x*+‘*,  x^  .  x"-*,  x"-*+*  (93) 

constitute  an  extended  Hermite  system  on  the  interval  (0, 1]  unth  a  an  arbitrary  non¬ 
integer  reed  number,  (see  Example  2.S  and  Definition  5.2) 


Theorem  5.1  Suppose  that  functions  {9i,9at' *  *  )92n}  Are  all  integrable  on  [a,  6],  and 
constitute  a  Chehyshev  system  on  (a,  6].  Suppose  further  that  on  the  interval  [a  +  6,5], 
the  n-point  Gaussian  quadrature  (7)  is  given  by  the  nodes 


(94) 

and  the  weights 

(95) 

for  anyS  €  (0,6  — a). 

Then  for  any  e>  0,  there  exists  Su>0  such  that  for  all  6  <  So, 

/  w(x)9,(x)dx-2u;fV<(»f^) 

<  e. 

(96) 

i=i 


for  oil  t  =  1,2,*  **  ,2n. 


f  u;(x)9<(x)dx  =  f;u>fV.(*5*^) 


Proof:  Due  to  the  Definition  2.4  of  the  Gaussian  Quadratures,  we  have 

r' 

/a+5 

for  all  t  =  1, 2,  *  *  ■ ,  n.  Subtracting 

J  w(x)9,(x)dx 

from  both  sides  of  (97),  we  obtain 

I  u{x)(pi{x)dx  -Y^wf^(pi{xf^)  =  [  u{x)(pi{x)dx  -  f  a;(x)9<(x)dx 
Ja  .■_.i  Ja  Ja+S 

/a+4 

u(x)(pi(x)dx. 


(97) 


(98) 
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Due  to  the  assumption  that  ipi  are  integrable  on  the  interval  [a,  fr]  for  all  ==  1, 2,  •  •  • ,  2n, 
for  any  e  >  0,  there  exists  ^  >  0  such  that  for  any  £  <  6o, 


IT' 


Lj{x)(pi(x)dx 


<  e 


(99) 


for  all  =  1,2, ••-,2n. 

Now,  (%)  follows  immediately  from  (98),  (99).  ■ 


Theorem  5.2  Suppose  that  under  the  assumptions  of  Theorem  2.2  and  Theorem  5.1, 

V>\  =  1-  (100) 

Then  there  exists  a  uni^e  n-point  Gaussian  quadrature  (7)  with  respect  to  the  functions 

»9an  (101) 

such  that  all  the  nodes  xi,X3, "•,Xn  lie  in  the  open  interval  (a, 6),  and  all  the  weights 
toi,  tPa,  •  •  • ,  u>n  ore  positive.  Furthermore,  for  alli=:  1, 2,  •  •  • , n, 

Umxj^  s=  Xi  (102) 

and 

=  Wi  (103) 

where  the  nodes  xj^^  and  the  weights  are  defined  in  Theorem  5.1. 

Proof:  Due  to  Theorem  2.2,  there  exists  a  uni<  ^ue  Gaussian  quadrature  (7)  with  respect 
to  the  functions  (101)  on  the  interval  [a, 6]  such  that  all  the  nodes  xi, xj,  •  •  • , x,  lie  in 
the  open  interval  (a,  6),  and  all  the  weights  tox,  tPa,  *  *  * ,  Wn  are  positive. 

On  the  other  hand,  due  to  Theorem  2.1,  for  any  ^  €  (0, 6  —  a),  there  exists  a  unique 
Gaussian  quadrature  (7)  with  respect  to  the  functions  (101)  on  the  interval  [a  +  5,  b] 
such  that  the  nodes  x|^^  €  (o  +  ^,  6),  and  the  weights 

>  0  (104) 

for  all  t  =  1, 2,  ”  * ,  n.  Combining  (100)  with  Definition  2.4  of  the  Gaussian  quadratures, 
we  obtain 

[  u;(x)v?i(x)dx  =  (105) 

Now,  for  any  ^  €  (0, 5  —  o),  we  will  define  two  vectors  xs,we  €  via  the  formulae 


where  and  u;j*^  are  the  Gaussian  nodes  and  weights  with  respect  to  functions  (101) 
on  the  interval  [a  +  6]. 

Clearly,  for  any  €  (0, 6  —  a),  we  have 


||x<l|  <c->/n 


with  c  3=  max(|a|,  |6|). 

Combining  (105)  and  (104),  we  obtain 

llwill  <d  y/n 

for  any  ^  €  (0, 6  —  a),  with  d  given  by  the  formula 


(108) 


(109) 

(110) 


We  will  show  that  for  all  £  €  (0,  b  —  a),  there  exists  only  one  limit  point  for  the  set 
of  vectors  x«,  and  only  one  limit  point  for  the  set  of  vectors  w«. 

Suppose  that  there  exists  a  sequence  of  positive  real  numbers  ^i,  ^,  *  *  * ,  such  that 
lim*_oo  Sk  =  0, 


Umxtf*  =  y, 

lim  wjj  =  V, 


(111) 
(112) 

with  y  =  (yi,ya,  •  •  •  fVnf,  and  v  =  (vi, vj, •  •  • , Wn)^. 

Due  to  Theorem  5.1,  Definition  2.4  of  the  Gaussian  Quadratures,  and  conditions 
(111),  (112),  we  have 


/  u{x)(pi{x)dx  =  lim  /  ui{x)(pi{x)dx 
Ja  k-*co  Ja+Si, 


=  Km;Cu;f*V.(#^) 

n 

=  Sw(yi) 


(113) 


for  all  t  =  1,2,- •  *,271. 

Due  to  (113)  and  the  uniqueness  of  the  Gaussian  quadrature  (7)  with  respect  to  the 
functions  (101)  on  the  interval  [a,  6],  we  have 

Vj  =  Xj  (114) 

Vj  =  Wj  (115) 


for  all  j  =  1,2,*' -,71,  where  xi,X2, ■**,Xn  are  the  nodes  and  101,103, ••*,tc;n  are  the 
weights  of  the  Gaussian  Quadrature  (7)  with  respect  to  the  functions  (101)  on  [a,  b].  In 
other  words,  the  set  of  vectors  x«  for  all  ^  €  (0,  b  —  a)  has  a  unique  limit  point 

X0  =  (Xi,X2,--*,Xn)^, 
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and  the  set  of  vectors  w«  for  all  ^  €  (0,  b  —  a)  has  a  unique  limit  point 

Now,  the  formulae  (102)  and  (103)  follow  from  (108),  (109),  and  the  fact  that  each 
of  the  two  sets  of  the  vectors  X(  and  w;  possesses  only  one  limit  point  in  i2".  ■ 

Remark  5.1  Clearly,  the  condition  =  1  tn  Theorem  5.2  can  he  relaxed.  To  insure 
the  boundedness  of  the  Gaussian  weights,  we  only  need  to  impose  the  condition  that  there 
exists  a  real  e  >  0,  and  2n  real  numbers  ai,  a^,  *  *  *,  oan  such  that 

^ai‘(pi{x)>e  (116) 

for  all  x  €  (a,  6). 

6  The  Numerical  Algorithm 

We  can  now  compute  Gaussian  quadratures  for  both  smooth  functions  and  functions 
with  end-point  singularities  using  the  numerical  apparatus  developed  in  Sections  4  and 
5.  The  Gaussian  quadrature  rules  for  an  extended  Hermite  system  can  be  obtained  by 
solving  a  system  of  non-linear  equations  (53).  Due  to  Theorems  4.5  and  5.2,  the  modified 
Newton’s  method  ddined  by  the  formula  (^)  converges  quadratically  when  it  is  applied 
to  the  system  of  equations  (53). 

As  is  well-known,  the  Newton  method  is  sensitive  to  the  choice  of  the  initial  approx¬ 
imation  Xo,  and  we  use  the  continuation  method  (see  Subsection  3.3  above)  to  obtain 
the  latter.  More  specifically,  given  an  extended  Hermite  system 

(117) 

on  the  interval  [a,  6],  we  construct  a  family  of  extmided  Hermite  systems 

(118) 

with  t  €  [0, 1],  and  such  that 

y»?(x)  =  (119) 

^J(x)  =  (pi{x),  (120) 

for  all  X  €  [o,  b],  and  i  =  1,2,  ••  •  .2n.  For  each  t  €  [0, 1],  we  construct  the  system  of 
equations  (53)  corresponding  to  the  extended  Hermite  system  (118)  via  the  formulae 
(39),  (40).  Clearly,  for  the  extended  Hennite  system  (119),  the  solution  of  the  system 
of  equations  (53)  is  known  (see  Remark  2.2),  and  we  use  the  continuation  method  (see 
Subsection  3.3)  to  obtain  the  solution  of  (53)  for  the  Hermite  system  (117). 

In  the  numerical  examples  of  the  following  section,  the  one-parameter  families  of 
Hermite  systems  are  constructed  as  follows. 
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1.  For  Hermite  systems  of  the  form 


(121) 

the  one-parameter  family  of  systems  is 

(p\{x)  =  (1  - 1)  •  x*~*  -I- 1  ■  (Pi{x)  (122) 

for  t  =  0, 1,*  •  •  ,2n  —  1. 

2.  For  Hermite  systems  of  the  form 

l,lnz,x,xlnx,x^,x^lnz,  •••  Inx,  (123) 

the  one-parameter  family  of  systems  is  given  by  (122). 

3.  For  Hermite  systems  of  the  form 

1, 5(x),  X,  xs(x),  x^  x^s(x),  •  •  • ,  x"”^,  x"“'s(x),  (124) 

with 

s(x)  =  x“,  (125) 

and  a  an  arbitrary  non-integer  real  ntunber,  the  one-paramet»  family  of  systems 
is 

1,  s{t,  x),  X,  X  •  s(t,  x),  x’,  x’  •  s(t,  x),  •  •  • ,  x""^ ,  x""^  •  s(t,  x),  (126) 

with 

s(t,x)  =  x*‘*.  (127) 

Remark  6.1  The  necessary  number  m  of  steps  in  the  continuation  process  (see  Subsec¬ 
tion  S.S)  is  significantly  reduced  if,  prior  to  the  application  of  the  above  procedure,  the 
original  system  (117)  is  orthonormalied  (for  example,  via  the  Gram-Schmidt  process). 
To  do  that,  we  discretize  the  original  functions  <pi  at  nested  Chebyshev  nodes  (see  Subsec¬ 
tion  S.  1 ),  and  perform  the  Gram-Schmidt  procedure  on  the  obtained  finite-dimensional 
representations. 

The  following  is  the  formal  description  of  the  numerical  algorithm  (excluding  the 
continuation  process). 

Initialization 

Comment  [  Build  the  structure  for  integration  and  interpolation.  ] 

Step  1 
do 

Subdivide  [a,  b]  into  subintervals  clustering  near  end-points  (see  Figure  3.1). 
enddo 
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Orthogonalization  (optional) 

Comment  [  Perfonn  Gram-Schmidt  orthogonalization  on  the  given  set  of  functions.  ] 

Step  2 


do  t  =  1,2, •••,2n 

do  j  =  1,2,  •••,*  —  1 

Orthc^onalize  the  t-th  function  (pi  with  respect  to  j-th  function  pj 

endo 

enddo 

Nested  Chel^shev  Apprcndmation  (optional) 

Comment  {  Generate  the  nested  Chebyshev  eq>ansions  for  the  orthogonalized  functions 
(see  Subsection  3.1).  ] 


Step  3 


do  t  =  1,2,*  •  ’  ,2n 

Construct  the  local  Chebyshev  expansion  of  the  t-th  function  <pi  based  on 
Observation  3.1. 
enddo 

Newton’s  Iteration 

Comment  [  Conduct  Newton’s  iteration  to  find  Gaussian  nodes  and  weights.  ] 

Step  4 


do 

Construct  functions  <r,'  and  t/i  for  t  =  1, 2,  • '  • ,  n  via  formulae  (39)  and  (40). 
enddo 


Step  5 


do 

Adjust  Gaussian  nodes  {x,}  via  formulae  (90). 
enddo 


Step  6 


do 

Compute  error  =  ||xfc+i  —  Xfc||. 
If  error  >  e,  Go  to  Step  4. 
enddo 
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Remark  6.2  The  procedure  described  above  requires  the  construction  of  the  functions  <Ti, 
qi,  i—  1, 2,  •  •  • ,  n,  given  the  functions  ipi,  i  —  1, 2,  •  •  • ,  2n.  The  latter  is  possible  for  any 
extended  Hermite  systems,  and  is  equivalent  to  inverting  the  matrix  (3).  Obviously,  for 
many  choices  of  functions  v’l,  v?2,  *  *  * ,  <fi2nt  the  matrix  (S)  will  be  ill-conditioned,  including 
the  numerical  examples  given  in  the  following  section.  Thus,  in  order  to  obtain  the  double 
precision  results  presented  in  this  paper,  the  authors  have  performed  all  computations  in 
extended  precision  (REAL  *Z2). 

7  Numerical  Results 

We  have  implemented  the  numerical  algorithm  described  in  Section  6  for  the  computa¬ 
tion  of  Gaussian  quadrature  formulae,  and  tested  it  on  various  example. 

Example  7.1  Gaussian  Quadratures  voith  respect  to  the  Bessel  Functions 

Jo,Ju“,J2n-\  (128) 

on  [0,10]  are  given  in  Table  1,  and  tested  on  selected  functions  in  Table  12. 

Example  7.2  Gaussian  Quadratures  with  respect  to  the  Bessel  Functions  (128)  on 
[0, 10]  with  the  weight  function 

(129) 

are  given  in  Table  2,  and  tested  on  selected  functions  in  Table  IS. 

Example  7.3  Gaussian  Quadratures  with  respect  to  the  system  of  functions 

l,lnx,a;,xlnx,x^,x^lnar,  •  •  •  ,x’*“^,x’‘~^  Inx  (130) 

on  [0, 1]  are  given  in  Table  3,  and  tested  on  selected  functions  in  Table  ll- 
Example  7.4  Gaussian  Quadratures  with  respect  to  the  systems  of  functions 

1,  x“,  X,  x'+“,  x^,  x2+“, . . . ,  x’*-^  x’’-^+“  (131) 

on  [0,1]  are  given  respectively  in  Tables  4-11  for 

2  1111112 
^  ~  3’  2’  3’  4’  4’  3’  2’  3’ 
and  tested  on  selected  functions  in  Tables  15-22  respectively. 

Remark  7.1  Systems  of  the  form  (128)  are  often  encountered  in  physics.  It  turns  out 
that  the  system  of  functions  (128)  on  the  interval  [0,5]  is  an  extended  Hermite  system 
only  for  certain  combinations  of  B  and  n.  A  somewhat  subtle  analysis  shows  that  the 
system  (128)  is  an  extended  Hermite  system  on  the  interval  [0, 5]  as  long  as  j  <  n  <  B. 
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8  Conclusions 


A  numeric&l  algorithm  has  been  presented  for  the  construction  of  the  generalized  Gaus¬ 
sian  quadrature  rules,  introduced  in  [6].  The  quadrature  rules  of  this  paper  possess  most 
of  the  desirable  properties  of  the  classical  Gaussian  integration  formulae,  such  as  p(^i- 
tivity  of  the  weights,  rapid  convergence,  mathematical  elegance,  etc.  The  algorithm  is 
applicable  to  a  wide  class  of  functions,  including  smooth  functions  (not  necessarily  poly¬ 
nomials),  as  well  as  functions  with  end-point  singularities,  such  as  those  encountered  in 
the  solution  of  integral  equations,  complex  analysis,  potential  theory,  and  several  other 
areas. 
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Table  1:  Gaussian  Quadrature  for  Bessel  Functions  of  the  First  Kind 


^10 

/  Jk-\{x)dx 
Jo 


N 


«=1 


for  fc  =  1,2,  •  •  ■  ,2N 


N 

Nodes  X, 

Weights  Wi 

5 

0.469238675868960E+00 

0.223157952970870E+01 

0.473407702933183E+01 

0.735478272434508E+01 

0.940238197915203E+01 

0.117179089779279E+01 

0.224136849121358E+01 

0.266214797121592E+01 

0.247001516625585E+01 

0.145415103193520E+01 

10 

0.1305356961 70244E+00 
0.672886387019932E+00 
0.159143208236292E+01 
0.280041847052746E+01 
0.419499640127942E+01 
0.566066246651666E+01 
0.707810341253441E+01 
0.8326218709542G4E+01 
0.928767980905348E+01 
0.986125239933237E+01 

0.333260223918652E+00 

0.742501752317741E+00 

0.108005279782693E+01 

0.132028812483455E+01 

0.144971022336955B+01 

0.146169692583314E+01 

0.135292016838366E+01 

0.112355472887245E+01 

0.782143625547308E+00 

0.353871429096005E+00 

25 


Table  2:  Gaussian  Quadrature  for  Bessel  Functions  with  Weight  Function 


/lO  I  N 

/  —7=Jk-iix)dx  -  2^  WiJk-iixi)  for  fc  =  1, 2,  •  •  ’ , 2JV 
Jo  yJX 


Weights  uji 


0.185798994260858E+01 

0.162483079218151E+01 

0.132302463248861E+01 

0.990971542854905E+00 

0.527439013034203&I-00 


5  0.221525014168133E+00 

0.181515943791217E+01 
0.43475721 1490782E+01 
0.710741692330460E+01 
0.933190129335545E+01 


0.586135127137856E-01 

0.517839232816169E+00 

0.138739079186343E+01 

0.258047692418516E+01 

0.398718764645530E+01 

0.548586592652336E+01 

0.694907018937755E+01 

0.824679794056492E+01 

0.925206018774378E+01 

0.985409749058602E+01 


0.966145917958189E+00 
0.939662007448660E+00 
0.889891823703254E+00 
0.821395494863014E+00 
0.737970249083253E+00 
0.641493179529387E+00 
0.531905946105306E+00 
0.408087704961759E+00 
0.269503595251694E+00 
0.1 18499400771494E+00 
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Table  3:  Gaussian  Quadrature  for  Products  of  Polynomials  and  Logarithmic  function 


N 

/  (pkix)dx  =  V Wi<pk{xi)  for  k  =  l,2,  --,2N 
where  {^,}  =  {l,lnar,x,xlnx, •••.x^“^,a:^'''lnx} 


N 

Nodes  X{ 

Weights  Wi 

5 

0.565222820508010E-02 

0.734303717426523E-01 

0.284957404462558E+00 

0.619482264084778E+00 

0.915758083004698E+00 

0.210469457918546E-01 

0.130705540744447E+00 

0.289702301671314E+00 

0.350220370120399E+00 

0.208324841671986E+00 

IT 

0.482961710689630E-03 

0.698862921431577E-02 

0.3261 13965946776E-01 

0.928257573891660E-01 

0.198327256895404E+00 

0.348880142979353E+00 

0.530440555787956E<fOO 

0.71676464851 1655E+00 

0.875234557506234E+00 

0.975245698684393E+00 

0.183340007378985E-02 

0.134531223459918E-01 

0.404971943169583E-01 

0.818223696589C36E-01 

0.129192342770138E+00 

0.169545319547259E+00 

0.189100216532996E+00 

0.177965753961471E4-00 

0.133724770615462E+00 

0.628655101 770325E-01 

15 

0.105784548458629E-03 

0.156624383616782E'02 

0.759521890320709E-02 

0.228310673939862E-01 

0.523886301568200E'01 

0.100758685201213E+00 

0.170740768849943E+00 

0.2625912061 18993E+00 

0.373536505184558E+00 

0.497746358414533E+00 

0.626789031392373E+00 

0.750516103461408E+00 

0.858255335207861E+00 

0.940141291212346E+00 

0.988401595986342E+00 

0.403217724648460E-03 

0.306297843478700E-02 

0.97842121 187661 5E-02 

0.215587522255813E-01 

0.383230673708892E-01 

0.588981990263004E-01 

0.811170299392595E-01 

0.102122101972069E+00 

0.118789059030401E+00 

0.128210316446694E+00 

0.128163327417093E+00 

0.117489465888492E+00 

0.963230185695904E-01 

0.661345398318934E-01 

0.296207140035355E-01 
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Table  3:  Gaussian  Quadrature  for  Products  of  Polynomials  and  Logarithmic  function 

,1  N 

/  (pkix)dx  =  ^i<Pk(xi)  for  A:  =  1, 2,  •  •  • ,  27V 

•'0  i=l 

where  {^i}  =  {l,ln*,x,arlnx,---,x^“^x^“*lnx} 


IP 

Nodes  Xi 

Weights  Wi 

"20“ 

0.352330453033401E-04 

0.52609398251 7410E-03 

0.258751954058141E-02 

0.793447194838041E-02 

0.186828881374457E-01 

0.370976733697505E.01 

0.653124886740214E-01 

0.10504850471 1551E+00 

0.157359691819002E+00 

0.222430062767455E+00 

0.299443765654100E+00 

0.386542446943882E+00 

0.480876453826790E+00 

0.578747932205507E+00 

0.675835475840038E+00 

0.767482460872564E+00 

0.849025253970320E+00 

0.916133703241664E+00 

0.965135427900256E+00 

0.993303536456954E+00 

0.134499676467758E-03 
0.103477692295062E-02 
0.337726367723322E-02 
0.767355619359468E-02 
0.142054962855420E-01 
0.229844384632086E-01 
0.337363605577136E-01 
0.459147630734522E-01 
0.587404799428040E-01 
0.71265013161 1020E-01 
0.824518089775832E-01 
0.912682015163873E-01 
0.967797159091613E-01 
0.982381433400897E-01 
0.951553030540297E-01 
0.873556504104574E-01 
0.750027772122717E-01 
0.585972958082337E-01 
0.3894725054961 14E-01 
0.171372052681059E-01 
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Table  3;  Gaussian  Quadrature  for  Products  of  Polynomials  and  Logarithmic  function 


/  (pk{x)dx  =  22  ^i<Pk{xi)  for  fc  =  1, 2,  •  •  • ,  2JV 
where  =  {l,in*,x,xlnx,- •  •  ,x^“^,x^“'lnx} 


N 

Nodes  Xi 

0.148805205646328E-04 

0.223091159576411E-03 

0.1 10464364905333E-02 

0.341946946887918E-02 

0.815052929502022E-02 

0.164289374947747E-01 

0.294459835598307E-01 

0.483575697078870E-01 

0.741870939196732E-01 

0.107732955883456E+00 

0.149486638258006E+(X) 

0.199566730959199E+00 

0.257673355831231E+00 

0.323066266406625E+00 

0.394568512069093E+00 

0.470596049553319EH-00 

0.549212146433099EH-00 

0.628203942243359E+00 

0.705177201069257E+00 

0.777664184415768E+00 

0.843238762138540E+00 

0.899S32416106158E+00 

0.9^M844733405702E+00 

0.977242575226688E+00 

0.995647215456440E+00 

0.568460660250201^04 

0.439997585768285E-03 
0.145071890475698E-02 
0.334401873816821E-02 
0.630809954735095E-02 
0.104488723103430E4)1 
0.157795036631243E-01 
0.222157908473636E-01 
0.295777024140889E>01 
0.375970456071 727E-01 
0.459308515949728E-01 
0.541797236656935E-01 
0.619100915223039E-01 
0.686790748928476E-01 
0.740604961651057E-01 
0.776705045127697E-01 
0.791912877674640E-01 
0.783914525088262E-01 
0.751418416138657E-01 
0.694258212157761E-01 
0.613433919048328E-01 
0.511088512980136E-01 
0.390421895640177E-01 
0.255554713626385E-01 
0.111503547267104E-01 

Table  3:  Gaussian  Quadrature  for  Products  of  Polynomials  and  Logarithmic  functi<»i 


/  (pk{x)dx  =  V  Wiifik{xi)  for  fc  =  1, 2,  •  •  • ,  2JV 
Jo  Srt 

where  {^j}  —  {l,ln«,x,xlna5,-*-,r^~*,x^“^lnx} 


Nodes  X,- 


0.732379743551900E-05 

0.110044700353982E-03 

0.5469183257031 79E-03 

0.170185751774368E-02 

0.408386360682336E-02 

0.830004117175941E-02 

0.150229781480879E-01 

0.249539236043865E-01 

0.387833861559789E-01 

0.571508984622219E-01 

0.806057414498694E-01 

0.109570394208134E4-00 

0.144308372971936E+00 

0.184897949395368E+00 

0.231213001514211E+00 

0.282911960162101E+00 

0.339435481155556E+00 

0.4000131 13074849E+00 

0.463678856906249E+00 

0.529295142710295E+00 

0.595584395297873E+00 

0.661 167040372616E+00 

0.724604528155535E+00 

0.78444573471 8383E+00 

0.839274951465986E+00 

0.887759597608297E+00 

0.928695795957382E+00 

0.9610500591841 75E+00 

0.98399570351 9948E+00 

0.996945958679506E+00 


0.279892154036191E-04 
0.217365526303388E-03 
0.720703585941237E-03 
0.167446096386701E-02 
0.319128240452640E-02 
0.535378831096153E-02 
0.820962136493378E-02 
0.117680292095091E-01 
0.159981435010663E-01 
0.208290410897408E-01 
0.261515976575626E-01 
0.318220682660372E-01 
0.376672559968915E-01 
0.434910622609641 E-01 
0.490821532269183E-01 
0.542224286606307E-01 
0.586959442912318E-01 
0.622979181160659E-01 
0.648434457091797E-01 
0.661755598539193E-01 
0.661722952804487E-01 
0.647524589265887E-01 
0.618798583538578E-01 
0.575658036673855E-01 
0.518697691962532E-01 
0.448981784990110E-01 
0.368013625032992E-01 
0.277688703796847E-01 
0.180238737856617E-01 
0.782767019615502E-02 


Table  4:  Gausnao  Quadrature  f(»r  Products  of  Polynomials  and  FVacti<mal  Powers 


/I  A 

I  <^k(x)<ix  =  Wi<pk(xi)  for  fc  =  1 , 2,  •  •  • ,  2iV 

''®  tal 

where  {^,}  =  {l,s(x),x,xs(x),---,x^“*,x^"^s(x)}  with  s(x)  =  xt 


Nodes  X,- 


0.11 1388121461 113E-01 
0.989954782999841E-01 
0.325462965706881E+00 
0.650376177175503E+00 
0.92383014138331 lE+00 


Weights  to,- 


0.350341916241438E-01 

0.152986023564027E+00 

0.293439234461264E+00 

0.329423482895757E+00 

0.189117067454808E+00 


10  0.106075936690850E-02 
0.105957835374351E-01 
0.419569285287569E-01 
0.109001951813403E+00 
0.219861071532861E+00 
0.372071233678686E+00 
0.550779272731362E+00 
0.730779332479715E+00 
0.881954366270197E+00 
0.976639980362546E+00 


0.342465634725548E-02 

0.179697406786380E-01 

0.471 712377689652E-01 

0.883592422887387E-01 

0.132980363859702E+00 

0.168929398195588E+00 

0.184191248381816E+00 

0.170661619172872E+00 

0.126951038006216E+00 

0.593614553002092E-01 


0.241818436310427E-03 

0.247904422657297E-02 

0.102231505109504E-01 

0.280798632789153E-01 

0.608339281242860E-01 

0.1 12443594490995E+00 

0.185105777705090E+00 

0.278545764961840E+00 

0.389652905586201E+00 

0.512532528698514E+00 

0.638982003700951E+00 

0.759334800695944E+00 

0.86356035521 1660E+00 

0.942468575671709E+00 

0.988861955722520E+00 


0.785216839155443E-03 

0.428667428787457E-02 

0.119847435540932E-01 

0.245432960884266E-01 

0.416404983681 140E-01 

0.619571219761760E-01 

0.833277026033653E-01 

0.103032984392016E+00 

0.118190780445325E+00 

0.126187019062355E+00 

0.125082141394901E+00 

0.113931212165169E+00 

0.929687707138682E-01 

0.63631 1469832558E-01 

0.28450691 1259040E-01 
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Table  4:  Gauanan  Quadrature  for  Products  of  Polynomials  and  FVactional  Powers 


/  (Pk(x)dx  =  5^  u)i(pk{xi)  for  fc  =  1, 2,  •  •  • ,  2iV 
where  with  s(x)  =  xi 


N 

Nodes  X,- 

Weii^ts  VDi 

IT 

0.822258720897910E-04 

0.851231 173182477E-03 

0.35641 1730804281E-02 

0.990446520349904E-02 

0.222308791404550E-01 

0.424339835363049E-01 

0.725741080670503E-01 

0.114174415866984E+00 

0.168083972212092E+00 

0.234302815001518E+00 

0.311876663295677E+00 

0.398872872080255E+00 

0.492442140014574E+00 

0.588962927640006E+00 

0.684258212417279E+00 

0.773867742847468E+00 

0.853353935948184E+00 

0.918616478895035E+00 

0.966190178403627E+00 

0.993508545383756E+00 

0.267576170194858E-03 

0.148254407742955E^)2 

0.424221866790033E-02 

0.897068317712191E-02 

0.158666804741713E-01 

0.248722570057384E-01 

0.356656880963368E>01 

0.476793965125403E4)1 

0.601415054233364E-01 

0.721375449972604E-01 

0.826870082942195E>01 

0.908280936995649E-01 

0.957031990037015E-01 

0.966376154468787E>01 

0.932044202436777E-01 

0.852697379326592E-01 

0.730142442020013E-01 

0.569289535722086E-01 

0.377866171625958E-01 

0.166140158404627E-01 
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Table  5:  Gaiunan  Quadrature  for  Products  of  PolyuOTuals  and  FVactional  Powers 


A 

I  (pk{x)dx  =  2^  toupkixi)  for  fc  s=  1, 2,  •  •  • ,  2iV 

ml 

where  {^i}  s=  {l,s(a5),*,xj(*),-*-,x^~*,x^“*5(x)}  with  «(x)  =  xi 


N 


5  0. 


laffiTl 


10 

0 

0 

0 

0 

0 

0 

0 

0 

0 


15  0 
0 
0. 
0 
0 
0. 
0. 
0 
0 
0 
0. 
0. 
0. 
0. 
0. 


Nodes  Xi 


970916313338209 
927420088040289E-01 
315872313916462E+00 
643182477910772E+00 
921965110615521E+00 


901742772555592E-03 

9660729921 18868E-02 

396093898716370E>01 

10501 1991918026E+00 

214610971 190650E+00 

366460914978464E+00 

545885024355929E+00 

727418879329945E+00 

880346704943949E+00 

976306802645093E+00 


203617338486320E-03 

223725122307619E>02 

954788886412147E-02 

267548306484035E-01 

587263787087472E-01 

109550941 121055E+00 

181570447192808E+00 

274636464058843E+00 

385717152459535E+00 

508930876164803E+00 

636017805747573E+00 

7571940440891 15E+00 

862273903642310E+00 

941904608637681E+00 

988750439427061E+00 


0.314958290433846E-01 

0.147817740145233E+00 

0.292773974169340E+00 

0.334349276188739E+00 

0.193563180453303E+00 


0.299828120481279E-02 

0.168386395659664E-01 

0.455491829065261E>01 

0.868038128143013E-01 

0.132106151 126701E+00 

0.169114219381655E+00 

0.185393787355447E+00 

0.172422600578352E+00 

0.1285743O9O18165E+OO 

0.601990160480740E-01 


0.680684768275793E-03 

0.397246938629008Er02 

0.114344447894424E-01 

0.238083807101952E-01 

0.408326185452397E-01 

0.612195391729371E-01 

0.828014703688452E-01 

0.102824724623773E+00 

0.1 183483841 19426E+00 

0.126687910323023E+00 

0.125836213297194E+00 

0.114797651013883E+00 

0.937832798881 199E-01 

0.642380897499318E-01 

0.287341392434250E-01 


Table  5:  Gaossiaii  Quadrature  for  Products  of  Polynomials  and  FVactional  Powm 


A 

/  (pk{x)dx  =  ^Wi<pk{xi)  for  = 

•'0  «i 

where  {^,}  =  with  3(x)  =  x^ 


N 

Nodes  Xt 

Wei^ts  to. 

w 

0.688907338392845E-04 

0.764138398652809E-03 

0.331021381556300E-02 

0.946809336603785E-02 

0.213344721137020EM)1 

0.410964086262251E-01 

0.707645876791935E-01 

0.111910522350377E+00 

0.165432926252467E+00 

0.2313761 14514019E+00 

0.308818842285795E+00 

0.395845863335002E+00 

0.489607131209166E+00 

0.586462066987943E+00 

0.682198172145692E+00 

0.772307249294583E+00 

0.852296694398426E+00 

0.918010354508321E+00 

0.965932757966444E+00 

0.993458519536807E+00 

0.230763814351417&03 
0.136612455217214E-02 
0.402291522792898E-02 
0.864643822538178E4)2 
0.154553358601049E-01 
0.244080922533623E-01 
0.351941337584140E4)1 
0.472505923470994£4)1 
0.598034836852591E-01 
0.719298383077259E-01 
0.826354136333313E-01 
0.909411778827062E-01 
0.959709232198866E-01 
0.970322101197066E-01 
0.936833837251 139E-01 
0.857805916759546E-01 
0.735004796732940E-01 
0.573364747582598E-01 
0.380699389274024E-01 
0.167416883525447E4)1 
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Table  6:  Gaussian  Quadrature  for  Products  of  Polynomials  and  Fractional  Powers 


[  (Pk{x)dx  =  53  ^ifpkixi)  for  fc  =  1 , 2,  •  • ,  2iV 
where  {(^,}  =  with  s(x)  =  xi 


N 

Nodes  X, 

Weights  Wi 

5 

0.831 174531456776E-02 

0.863966362795308E-01 

0.305943516943443E+00 

0.635656558720652E+00 

0.920004207024857E+00 

0.279778782123048E-01 

0.142398935990482E+00 

0.291943668689807E+00 

0.339447240627363E+00 

0.198232276480043E+00 

10 

0.751897878625465E-03 

0.874666094268371E-02 

0.372684363855664E-01 

0.100986759928202E+00 

0.209274931334610E+00 

0.360730157291068E+00 

0.540868549812902E+00 

0.723966527445448E+00 

0.878692655294032E+00 

0.975963746745065E+00 

0.259000938099737E-02 

0.157085175605557E-01 

0.438972012175236E-01 

0.851979518885884E-01 

0.131185462315937E+00 

0.169279428891491E+00 

0.186612247123971E+00 

0.174225330076043E+00 

0.130242598410538E+00 

0.610612531343552E-01 

15 

0.168121349224655E-03 

0.200427399422478E-02 

0.888445820257933E-02 

0.254380748331252E-01 

0.566163366479501E-01 

0.106639742234206E+00 

0.1 77998880856086E+00 

0.270675783247184E+00 

0.381720959785030E+00 

0.505267785653412E+00 

0.632999186646096E+00 

0.75501 1893747991E+00 

0.860961671393287E+00 

0.9413290761 90074E+00 

0.988636608057791E+00 

0.582014429388243E-03 

0.366362217710688E-02 

0.108842247988913E-01 

0.230661756216959E-01 

0.400108265603416E-01 

0.604643908138577E-01 

0.822581407444590E-01 

0.102603879506980E+00 

0.118500791 181315E+00 

0.127191992481750E+00 

0.126600784497532E+00 

0.1 15679059680403E+00 

0.946133646757434E-01 

0.648572778769515E-01 

0.290234549535849E-01 

Table  6:  Gaussian  Quadrature  for  Products  of  Polynomiab  and  Fractional  Powers 


/  <pk{x)dx  =  ^  Wi<pk{xi)  for  Ar  =  1, 2,  •  •  • ,  2JV 

ml 

where  {v>.}  =  {l,s(x),x,xs(*),---,ar^”*,x^“^s(x)}  with  s(x)  =  xi 


N 

Weights  Wi 

IcT 

Q.565910606963981E-04 

0.680854395380959E-03 

0.3062701 10956997E-02 

0.894917200513315E-02 

0.204442383726108E-01 

0.39761 1592475325E-01 

0.689512869857525E-01 

0.109635234283660E+00 

0.162762389006638E+00 

0.228422459407150E+00 

0,305728258549570E+00 

0.392782689780531E+00 

0.486735430814252E+00 

0.583926798874029E+00 

0.680108431067843E+00 

0.770723438587922E+00 

0.851223225892569E+00 

0.91 7394745239490E+00 

0.965671257181204E+00 

0.993407695194012E+00 

0.196275471368437E-03 
0.125263851254803E-02 
0.380561564887333E-02 
0.832217213031223E-02 
0.150414196754775E-01 
0.239388575134562E-01 
0.347155654274889E-01 
0.468137478429227E-01 
0.594574979240510E-01 
0.717153419416687E-01 
0.825791227715473E-01 
0.910522631915573E-01 
0.962395518442346E-01 
0.974303994336976E-01 
0.941680397096348E-01 
0.862983560861 lOOE-01 
0.739938065555806E-01 
0.577502214287376E-01 
0.3835771 14718547E-01 
0.1687139541 88783E-01 
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Table  7:  Gatissian  Quadrature  for  Products  of  Polynomials  and  Fractional  Powers 


fi  ^ 

/  (pkix)dx  =  22  ^i<fhixi)  for  fc  =  1, 2,  •  •  • ,  27V 

i=i 

where  =  {I,s(i),z,x5(a5), with  s(z)  = 


N 

Nodes  Xi 

Weights  to, 

~ 

0.762794972812507E-02 
0.831895488095584E-01 
0.3008440374531 26E+00 
0.631759839813656E+00 
0.9189851 15815450E+00 

0.26229803590671 7E-01 

0.139587984773339E+00 

0.291460348074849E+00 

0.342065172592006E+00 

0.200656690969134E+00 

10 

0.680661832514276E-03 

0.829805126031941E-02 

0.361008488950617E-01 

0.989606158324035E-01 

0.206573413000166E+00 

0.357817634747240E+00 

0.538312396928998E+00 

0.722204307474990E+00 

0.87784743501 1894E+00 

0.975788344331792E+00 

0.239309679083200E-02 

0.151439816123742E-01 

0.430595363165574E-01 

0.843751810440922E-01 

0.130706769333825E+00 

0.169354243901417E+00 

0.187227677434183E+00 

0.175143061430482E+00 

0.131094401625182E+00 

0.6150205051 10549E-01 

Is" 

0.151429644289982E-03 

0.189120816621145E-02 

0.855735146816403E-02 

0.247829453316383E-01 

0.555603949382247E-01 

0.105176990238943E+00 

0.176199068882164E+00 

0.268675523469233EH-00 

0.37969941 7736514E+00 

0.503412401869408E+00 

0.631468760396525E+00 

0.753904757831583E+00 

0.860295553378594E+00 

0.941036823465374E+00 

0.988578794355859E+00 

0.534953080073061E-03 

0.351 128622129104E-02 

0.106091654389163Er01 

0.226922601463664E-01 

0.395945355599134E-01 

0.600799922753622E-01 

0.819798274477196E-01 

0.102488560198090E+00 

0.118574970046103E+00 

0.127445268951420E+00 

0.126987144614300E+00 

0.116125573735615E+00 

0.950344508261 170E-01 

0.651716225755041E-01 

0.291703888832093E-01 

37 


Table  7:  Gaussian  Quadrature  for  Products  of  Polynomials  and  Fractional  Powers 


N 

I  <Pk{x)dx  =  ^Wi(pk{xi)  for  fc  =  l,2,---,2Ar 

i=i 

whCTe  {^i}  =  {l,s(x),x,xs(a:),*--,a:^~^,x^“^s(x)}  with  5(x)  =  xi 


N 

Nodes  Xt 

Weights  Wi 

20 

0.508401 144727475E-04 

0.640671409304506E-03 

0.294139627824379E-02 

0.869257229367654E-02 

0.200014921598775E-01 

0,390944230688757E-01 

0.680431707412333E-01 

0.108493186848796E+00 

0.161419592745879E+00 

0.226935224532033E+00 

0.304170319054039E+00 

0.391237157031920E+00 

0.485285420936866E+00 

0.582645888835788E+00 

0.679052102488331E+00 

0.769922539292854E+00 

0.850680233407198E+00 

0.917083282994612E+00 

0.965538933432666E+00 

0.993381975142296E+00 

0.179921229090492E-03 

0.119702077067209E-02 

0.369773784996860E-02 

0.816003283674897E-02 

0.148334720914000E-01 

0.237022859923791E-01 

0.344735758S59312E-01 

0.465922232807169E-01 

0.592814319750155E-01 

0.716054742949001E-01 

0.825491665123296E-01 

0.911070366015320E-01 

0.963742175415786E-01 

0.976308832970205E-01 

0.944125656226918E-01 

0.865599052197165E-01 

0.742432057258809E-01 

0.579594958287601E>-01 

0.385033140388101E-01 

0.169370334348568E-01 
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Table  8:  Gaussian  Quadrature  for  Products  of  Polynomials  and  Fractionad  Powers 


fi  " 

ipk{x)dx  =  ^wnpk{xi)  for  *:  =  l,2,-  -,2iV 

i=i 

where  {y?,}  =  {l,s(x),x,xs(x), •  •  with  s(x)  =  x~i 


N 

Nodes  Xi 

Weights  Wi 

5 

0.382423605041850E-02 
0.634597428923166E-01 
0.268079430001204E+00 
0.606207314918122E+00 
0.9122431 19495094E+00 

0.159917283936684E-01 
0.121052296479220E+00 
0.287406177837817E+00 
0.358886538601 875E+00 
0.216663258687420E+00 

10 

0.311664421263886E-03 

0.573840166641843E-02 

0.291436456086289E-01 

0.8659896401 95454E'01 

0.189850316327027E+00 

0.339616983180453E+00 

0.522237733374721E+00 

0.711076275414628E+00 

0.872496278362019E+00 

0.974676373102776E+00 

0.132408316145357E-02 

0.117668861542949E-01 

0.378545434895510E-01 

0.791317874471543E-01 

0.127550269062920E+00 

0.16%82299097476E+00 

0.191015829563274E+00 

0.180901690009248E+00 

0.136476978227289E+00 

0.642956337873400E-01 

15 

0.6715371 72803988E-04 

0.126348368805257E-02 

0.666285677053233E-02 

0.2089963924471 29E-01 

0.492092869824728E-01 

0.962910504069878E-01 

0.165187187294757E+00 

0.256372388322078E+00 

0.367215812312105E+00 

0.491920272820932E+00 

0.621967797979428E+00 

0.747019986367819E+00 

0.856148272384472E+00 

0.9392157995951  lOE+00 

0.988218403767777E+00 

0.286236608068052E-03 

0.262818581079834E-02 

0.895931713383214E-02 

0.204062530697822E-01 

0.370154044724167E-01 

0.576707391266183E-01 

0.802098718885372E-01 

0.101723142015267E+00 

0.118989956384134E+00 

0.128983991417795E+00 

0.129367050472822E+00 

0. 1 1 8892374347340E+00 

0.976520449023468E-01 

0.671291997815284E-01 

0.300862325687140E-01 
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Table  8:  Gaussian  Quadrature  for  Products  of  Polynomials  and  Fractional  Powers 


■1  N 

(f>k{x)dx  =  Wiipkixi)  for  fc  s=  1, 2,  • '  • ,  2JV 

I  i _ « 


where  {(^,}  =  with  s(x)  =  x“i 


N 

Nodes  X,- 

Weights  Wi 

0.221820235538413E-04 

0.4206681563881 07E-03 

0.224897834059227E-02 

0.719410590426075E-02 

0.173785665899570E-01 

0.351053065795165E-01 

0.625706338467844E<01 

0.101573342538239E+00 

0.153248856878474E+00 

0.217855193592975E+00 

0.294632990834409E+00 

0.381755192469280E+00 

0.47637381 7974157E+00 

0.574762338493956E+00 

0.672543338723645E+00 

0.764983189265488E+00 

0.847329128952469E+00 

0.915160094733383E+00 

0.964721591676480E+00 

0.99322307681 8039E+00 

0.946625535057872E-04 

0.879580435485498E-03 

0.306154371288216E-02 

0.718677600572896E-02 

0.135708222783955E-01 

0.222536266859334E-01 

0.329813226328586E-01 

0.452169532360338E-01 

0.581794488048750E-01 

0.709075209058289E-01 

0.823427811731510E-01 

0.914245608123939E-01 

0.971877635212259E-01 

0.988547376185605E-01 

0.959126372029696E-01 

0.881690173465919E-01 

0.757803722977461E-01 

0.592508853437043E-01 

0.394024605588481E-01 

0.173425268732812E-01 
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Table  9:  Gaussian  Quadrature  for  Products  of  Polyuomials  and  FVactional  Powers 


/■»  A 

/  ipkix)dx  =  ^  Wi(pk{xi)  for  A:  =  1 , 2,  ■  •  • ,  2iV 
.=1 

where  {y?,}  =  {l,s(x),x,xs(x), with  s(x)  =  x"^ 


N 

Nodes  Xi 

Weights  Wi 

5 

0.32569488S051994E-02 

0.600861836255002E'01 

0.262200689766926E+00 

0.601527262010120E-f00 

0.910998040648555E+00 

0.143436825038591E-01 

0.117635462047981E+00 

0.286500581347155E+00 

0.361906424437754E+00 

0.219613849663250E+00 

10 

0.261010670878138E-03 

0.533570406879994E-02 

0.279923598777842E-01 

0.8450010671 79312E-01 

0.186967540598809E+00 

0.336449186747933E+00 

0.519422386872603E+00 

0.709119391752190E+00 

0.871552940481438E+00 

0.974480099088773E+00 

0.116634009622830E-02 

0.112056554482295E-01 

0.369538395786576E-01 

0.782012224225048E-01 

0.126971878207205E+00 

0.169715037099809E+00 

0.191665171997203E+00 

0. 18190813623421 7E+00 

0.137424134774237E+00 

0.647885841417101E-01 

Is" 

0.559212686135403E-04 

0.116764741757783E-02 

0.635876913141328E-02 

0.202601345786029E-01 

0.481469972090806E-01 

0.947893055727041E-01 

0.163312504831184E+00 

0.254266678317515E+00 

0.365070714120304E+00 

0.489939651507906E+00 

0.620326692664507E+00 

0.745828824368843E+00 

0.855429899278146E+00 

0.938900129707038E+00 

0.988155905432853E+00 

0.250636999997100E-03 

0.248629682453753E-02 

0.868418049215750E-02 

0.200174031540757E-01 

0.365707842279019E-01 

0.572506624316614E-01 

0.798969454587593E-01 

0.101582504459800E+00 

0.119053915578673E+00 

0.129244090759107E+00 

0.129774994114533E+00 

0.119369417740098E+00 

0.981047682668204E-01 

0.674683579975495E-01 

0.302450414943277E-01 
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Table  9:  Gaussian  Quadrature  for  Products  of  Polynomials  and  Fractional  Powers 


f  (pk{x)dx  =  53  ^iV>hixi)  for  A:  =  1 , 2,  •  ■  • ,  2iV 

»=i 

where  {(^j}  =  {l,s(x),x,is(x),- •  •  with  s(x)  =  x“i 


N 

Nodes  X,- 

Weights  to, 

loT 

0.1841949075981 79E-04 
0.387587458269692E-03 
0.213955259577947E-02 
0.695120073310895E-02 
0.169467664404234E-01 
0.34441 74520001 76E-01 
0.616535209242469E-01 

0. 100407215362434E+00 
0.151866019233245E+00 
0.216313281828314E+00 
0.293009074582624E+00 
0.380137222473773E+00 
0.474850532149502E+00 
0.573412898781505E+00 
0.671427980742799E+00 
0.764136026657879E+00 
0.846753983774620E+00 
0.914829857705509E+00 
0.964581196697193E+00 
0.993195777887590E+00 

0.826457102725699E-04 
0.829421040741860E-03 
0.295732155968210E-02 
0.702434283841401E-02 
0.133575962599994E-01 
0.220068985356573E-01 
0.327254003368527E-01 
0.449795404557805E-01 
0.579877214786053E-01 
0.7078438381 19479E-01 
0.823037488978054E-01 
0.914756213036679E-01 
0.973244581548574E-01 
0.990625442531322E-01 
0.961685890414821E-01 
0.884443478050775E-01 
0.760438602249420E-01 
0.594724986507429E-01 
0.395568714306716E-01 
0.174121882096675E-01 
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Table  10:  Gaussian  Quadrature  for  Products  of  Polynomials  and  Fractional  Powm 


[  <fik{x)dx  =  Yi  ^i¥>k(xi)  for  fc  =  1 , 2,  •  •  • ,  2JV 

’'®  «=i 

where  with  s(x)  =  x~i 


N 

Nodes  Xj 

Weights  Wi 

5 

0.220055532702321E-02 
0.532526444285811E-01 
0.250000000000000E+00 
0.591 721954534264E+00 
0.908380401265687E+00 

0.111142584286221E-01 

0.110450910249386E+00 

0.284444444444444E+00 

0.368177760249980E+00 

0.225812626627567E+00 

10 

0.170217313506295E-03 

0.455197375232787E-02 

0.256945562245545E-01 

0.802601948484878E-01 

0.181 103722710989E+00 

0.329978061692620E+00 

0.513655588977735E+00 

0.705104124523579E+00 

0.869615340441312E+00 

0.974076745830678E+00 

0.869843410720295E-03 

0.100832309490842E-01 

0.351 184957693976E-01 

0.762838816843756E-01 

0.125764125553643E+00 

0.169760099161 llOE+00 

0.192982837625621E+00 

0.183967866746584E+00 

0.1393681 18201 496E+00 

0.658015008979678E-01 

15 

0.360449058720925E-04 

0.983656825228959E>03 

0.576031032998391E-02 

0.189863966971689E-01 

0.460162191690594E-01 

0.917631475619828E-01 

0.159522718866145E+00 

0.250000000000000E+00 

0.360716812863579E+00 

0.485914494639546E+00 

0.616988391777598E+00 

0.743404128057339E+00 

0.85396689374041  lE+00 

0.938257049225935E4-00 

0.988028562926358E+00 

0.184634499540016E-03 

0.220691 172454993E-02 

0.813303209689366E-02 

0.192316020292443E-01 

0.356670580668843E-01 

0.563926955430610E-01 

0.792541212080791E-01 

0.101289120962781E+00 

0.119177364119033E+00 

0.129768304472501E+00 

0.130602147750110E+00 

0.120339075896910E+00 

0.990261883702783E-01 

0.681591357635582E-01 

0.305686074965773E-01 
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Table  10:  Gaussian  Quadrature  for  Products  of  Polynomials  and  Fractional  Powers 


ri  ^ 

/  (pk{x)dx  =  ^  xoiifikixi)  for  fc  =  1, 2,  •  •  • ,  2JV 

Jo 

where  {v?«}  =  with  s(x)  =  x“i 


N 

Nodes  X, 

Weights  xvi 

~w 

0.118040372897702E-04 

0.324505506016988E-03 

0.192569889609293E-02 

0.647083718891321E-02 

0.160868754200355E-01 

0.331142308281794E-01 

0.598127724451431E-01 

0.980610158280346E-01 

0.149078672924258E+00 

0.213200816542450E+00 

0.289727337675948E+00 

0.37686452406d904E+00 

0.471767104543454E+00 

0.570679774395970E+00 

0.66916791 1554694E+00 

0.762418781880186E+00 

0.84558780901 1132E+00 

0.914160127147419E+00 

0.964296432783931E+00 

0.993140403222385E+00 

0.605164515048587E-04 
0.731395632734516E-03 
0.275022407735182E-02 
0.669890718081941E-02 
0.129282095841639E-01 
0.215082324328235E-01 
0.322066292668297E-01 
0.444969640416654E-01 
0.575967453908003E-01 
0.7053185091 11265E-01 
0.822215362195993E-01 
0.915762410818034E-01 
0.975991452767166E-01 
0.994820091823469E-01 
0.966862995286949E-01 
0.890019102330765E-01 
0.765778343958853E-01 
0.599218242567572E-01 
0.398700341676524E-01 
0.175534906876473E-01 
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Table  11:  Gaussian  Quadrature  for  Products  of  Polynomials  and  fractional  Powers 


/  (pk{x)dx  =  2^  Wi(pk(xi)  for  ir  =  1 , 2,  •  •  • ,  2JV 

.*1 

where  {^i}  =  {l,s(x),x,xs(*), with  s(x)  =  x~i 


N 

Nodes  Xi 

Weights  to, 

~ 

0.127054140407448E'02 

0.462681413226638E-01 

0.237083298080667E+00 

0.581208702007414E+00 

0.905561736397846E+00 

0.798095300782248E-02 
0.102702832369157E+00 
0.2820090891 12635E+00 
0.374825379887486E+00 
0.232481745622901E+00 

'W 

0.945849186124229E'04 

0.379603281250484E-02 

0.233952614888275E'01 

0,759462672091217E-01 

0.175082164491838E+00 

0.323295565014521E+00 

0.507675337823808E+00 

0.700934090961697E+00 

0.867600470065010E+00 

0.973657038209960E+00 

0.599207947722056E-03 

0.895661938070323E-02 

0.332290268550819E-01 

0.742813877203913E-01 

0.124481407759041E-H)0 

0.169777375189087E+00 

0.194331900792034E+00 

0.186099994751647E+00 

0.141387747529539E+00 

0.668553320747537E-01 

0.197855962731515E'04 

0.809581023980278E'03 

0.517321860359321E'02 

0.177160194789544E-01 

0.438708712542835E'01 

0.886977166657579E-01 

0. 155667671624679E+00 
0.245646859906903E+00 
0.3562649725S3447E+00 
0.481 792161982642E+00 
0.613565444840732E+00 
0.740915815910203E+00 
0.852464602157936E+00 
0.937596443363357E+00 
0.987897722986291E+00 

0.125556576654582E-03 

0.193257660045135E-02 

0.757902277249522E-02 

0.184325783933937E-01 

0.34741241 1366044&01 

0.555083582856597E-01 

0.785867325255755E-01 

0.100978745423403E+00 

0.119295279382622E+00 

0.130299495385019E+00 

0.131446824433812E+00 

0.121332382739473E+00 

0.9997161 76822775E-01 

0.688685445138014E-01 

0.309010441487584E-01 
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Table  11:  Gaussian  Quadrature  for  Products  of  Polynomials  and  FVactional  Powers 


/I  N 

/  <pkix)dx  =  22  for  A:  =  1 , 2,  •  •  • ,  2JV 

ial 

where  {^,}  =  {l,s(x),x,*s(x),- •  with  s(x)  *  x“t 


N 

Nodes  Xi 

Weights  Wi 

“20" 

0.644080669471020E-05 

0.265385497209493E-03 

0.171795135681289E-02 

0.599650565966228E-02 

0.152297986688364E-01 

0.317830016969808E-01 

0.5795901 15876318E-01 

0.956908382982334E-01 

0.146256150019543E+00 

0.210043256963358E+00 

0.286393214270049E+00 

0.373535741616497E+00 

0.468627945387644E+00 

0.567895206083029E+00 

0.666863963449413E+00 

0.760667402163021E+00 

0.844398042013510E+00 

0.913476675531666E+00 

0.964005785910473E+00 

0.993083879404685E+00 

0.408973476604909E-04 
0.636150271 129477E-03 
0.254437715881698E-02 
0.637191713688041E-02 
0.1249391471 70132E-01 
0.210015095348042E4)1 
0.316775070055125E>01 
0.44003051 1092655E-01 
0.571949651553830E-01 
0.702705143585106E-01 
0.821337066503373E-01 
0.916751595069523E-01 
0.978762191044640E-01 
0.999075512797006E-01 
0.972128770925840E-01 
0.895698602397595E-01 
0.771222579836698E-01 
0.603802128097188E-01 
0.401896290968344E-01 
0.176977224410030E-01 
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Table  12;  Integration  by  the  Gaussian  Quadrature  of  Bessel  Functions 

in  TABLE  1 

f(x)dx  »  53to,/(*.) 

ml 


/(x)  =  8in(x) 


/o**“/(x)dx  =  0.183907152907645EH 

>01 

N 

Computed  Integral 

Absolute  Error 

Relative  Error 

5 

0.183911205770184E+01 

0.40529E-04 

■rii>y7r!5T:rgiy 

10 

0.183907152907645E+01 

0.15543E-14 

0.84516E>15 

/(x)  =  -C06(x) 


/o"'/(x)dx  =  0.544021110889370E+00 

N 

Computed  Integral 

Absolute  Error 

Relative  Error 

5 

0.544074430126847E+00 

0.53319&04 

0.98010E-04 

10 

0.544021 1 10889362E+00 

0.77716E-14 

0.14285E-13 

fix)  =  -X 


Ji*" /(*)<!*  =  0.500000000000000EH 

N 

Computed  Integral 

Absolute  Error 

Relative  Error 

5 

0.499933635432992E+02 

0.66365E-02 

0.13273E-03 

10 

0.499999999999999E+02 

0.14921E-12 

0.29843E-14 
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Table  13:  Integration  by  the  Gaussian  Quadrature  of  Bessel  Functions 

in  TABLE  2 


/■*°  1  JL 

I  ^  £  WifiXi) 

•'®  V*  Si 


/(x)  =  sin(x) 


/o*";^/(x)dx  =  0.152512353028332E+01 

JV 

Computed  Integral 

Absolute  Error 

Relative  Error 

5 

0.152515304377709E+01 

0.29513E-04 

0.19352E-04 

10 

0.152512353028335E+01 

0.33307E>-13 

0.21839E’13 

/(*)  = 

1  —  cos(x) 

/o“;^/(x)dx  =  0.522924912143076E+01 

iV 

Computed  Integral 

Absolute  Error 

Relative  Error 

5 

0.522898009347896E+01 

0.26903E-03 

0.51447E-04 

10 

0.522924912142827E+01 

0.24931E-11 

0.47676E-12 

fix)  =  X 


/o*“;^/(x)dx  =  0.210818510677892E+02 

N 

Computed  Integral 

Absolute  Error 

Relative  Error 

5 

0.210781199174431E+02 

0.37312E-02 

0.17698E-03 

10 

0.210818510677890E+02 

0.19895E-12 

0.94371E-14 
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Table  14:  Integration  by  the  Gaussian  Quadrature  for  the  Products 
of  Polynomials  and  Logarithmic  Functions 
in  TABLE  3 


/  fix)dx  si'^wJixi) 


/(i)  =  sin(17x) 


f^  f{x)dx  =  0.750096081206822E-01 

N 

Computed  Integral 

Absolute  Error 

Relative  Error 

5 

-0.445245704584552E+00 

0.52026E+00 

0.69358E-1-01 

10 

0.759657791361190E-01 

0.9561 7E-03 

0.12747E-01 

15 

0.750095713818829E-01 

0.36739E-07 

0.48979E-06 

20 

0.75009608120871 lE-01 

0.18893E-12 

0.25188E-11 

25 

0.750096081206832E-01 

0.10131E-14 

0.13506E-13 

30 

0.750096081206823E-01 

0.11102E-15 

0.14801E-14 

fix)  = 


J^nx)dx  =  0,21 

i7391304347826E-01 

N 

Computed  Integral 

Absolute  Error 

Relative  Error 

5 

0.397077526634515E-02 

0.17768E-01 

0.81734E+00 

10 

0.20681633541 7969E-01 

0.10575E-02 

0.48645E-01 

15 

0.217347730977134E-01 

0.43573E-05 

0.20044E-03 

20 

0.217391290710919E-01 

0.13637E-08 

0.62730E-07 

25 

0.217391304347564E-01 

0.26222E-13 

0.12062E-11 

30 

0.217391304347826E-01 

0.69389E-17 

0.31919E-15 

fix)  =  — x^^lnx 


/o"  /(x)dx  =  0.472589792060492E-03 

N 

Computed  Integral 

Absolute  Error 

Relative  Error 

5 

0.349440397064675E-03 

0.12315E-03 

0.26058E+00 

10 

0.5543969482101 19E-03 

0.81807E-04 

0.17310E+00 

15 

0.473519740800504E-03 

0.92995E-06 

0.19678E-02 

20 

0.47259036015371 OE-03 

0.56809E-09 

0.12021E-05 

25 

0.472589792079166E-03 

0.18675E-13 

0.39515E-10 

30 

0.472589792060492E-03 

0.37947E-18 

0.80296E-15 
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Table  15:  Integration  by  the  Gaussian  Quadrature  for  the  Products  of 
Polynomials  and  Fractional  Power  s(x)  =  xi 
in  TABLE  4 


t=i 


/(x)  =  sin(lOx) 


fo*  fix)dx  =  0.183907152907645E+00 

N 

Computed  Integral 

Absolute  Error 

Relative  Error 

5 

0.205829062722746E+00 

0.21922E-01 

0.11920E+00 

10 

0.183906457937759E+00 

0.69497E-06 

0.37789E-05 

15 

0.183907152908513E+00 

0.86728E-12 

0.47158E-11 

20 

0.183907152907645E+00 

0.19429E-15 

0.10565E-14 

fix)  = 


fd  f(x)dx  =  0.384615384615385E-01 

N 

Computed  Integral 

Absolute  Error 

Relative  Error 

5 

0.261006235736059E-01 

0.12361E-01 

0.32138E+00 

10 

0.384350080122515E-01 

0.26530E-04 

0.68979E-03 

15 

0.384615380725112E-01 

0.38903E-09 

0.10115E-07 

20 

0.384615384615385E-01 

0.69389E-16 

0.18041E-14 

fix)  =  x^s{x) 


/o'  /(x)dx  =  0.377358490566038E.01 

N 

Computed  Integral 

Absolute  Error 

Relative  Error 

5 

0.261006235736059E-01 

0.12361E-01 

0.32138E+00 

10 

0.384350080122515E-01 

0.26530E-04 

0.68979E-03 

15 

0.3846153807251 12E-01 

0.38903E-09 

0.10115E-07 

20 

0.384615384615385E-01 

0.69389E-16 

0.18041E-14 
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Table  16:  Integration  by  the  Gaussian  Quadrature  for  the  Products  of 
Polynomials  and  Fractional  Power  s(x)  = 
in  TABLE  5 


f  f{x)dx  =  '^Wif{xi) 

«=i 


/(z)  =  sin(lOz) 


fi  fix)dx  =  0.183907152907645E+00 

N 

Computed  Integral 

Absolute  Error 

Relative  Error 

5 

0.205245178663791E+00 

0.21338E-01 

0.11603E+00 

10 

0.183906379065431E+00 

0.77384E-06 

0.42078E-05 

15 

0.183907152908629E+00 

0.98405E-12 

0.53508E-11 

20 

0.183907152907645E+00 

0.27756E-16 

0.15092E-15 

fix)  =  z»* 


/o"  fix)dx  =  0.384615384615385E-01 

N 

Computed  Integral 

Absolute  Error 

Relative  Error 

5 

0.253%7007271653E-01 

0.13065E-01 

0.33969E+00 

10 

0.38431 1761588644E-01 

0.30362E-04 

0.78942E-03 

15 

0.384615379738251E-01 

0.48771E-09 

0.12681E.07 

20 

0.384615384615383E-01 

0.20817E-15 

0.54123E-14 

fix)  =  z“s(z) 


/o"  fix)dx  =  0.377358490566038E-01 

N 

Computed  Integral 

Absolute  Error 

Relative  Error 

5 

0.243848070841558E-01 

0.13351E-01 

0.35380E+00 

10 

0.377002264605833E-01 

0.35623E-04 

0.94400E-03 

15 

0.3773584831 10791E-01 

0.74552E-09 

0.19756E-07 

20 

0.377358490566035E-01 

0.22898E-15 

0.60681E-14 
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Table  17:  Integration  by  the  Gaussian  Quadrature  for  the  Products  of 
Polynomials  and  Fractional  Power  3(x)  = 
in  TABLE  6 


[  f{x)dx  = 

i=i 


f{x)  =  sin(lOx) 


/o*  f(x)dx  =  0.183907152907645E+00 

N 

Computed  Integral 

Absolute  Error 

Relative  Error 

5 

0.203627012649713E+00 

0.19720E-01 

0.10723E+00 

10 

0.183906319056362E+00 

0.83385E-06 

0.45341E-05 

15 

0.183907152908719E+00 

0.10739E-11 

0.58396E-11 

20 

0.183907152907644E+00 

0.86042E-i5 

0.46786E-14 

fix)  = 


/o^  /(x)dx  =  0.384615384615385E-01 

N 

Computed  Integral 

Absolute  Error 

Relative  Error 

5 

0.246599205561 149E-01 

0.13802E-01 

0.35884E+00 

10 

0.384268064253340E-01 

0.34732E-04 

0.90303E-03 

15 

0.384615378508424E-01 

0.61070E-09 

0.15878E-07 

20 

0.384615384615384E-01 

0.11102E-15 

0.28866E-14 

/(x)  =  x*®s(x) 


si  f{x)dx  =  0.380228136882129E-01 

N 

Computed  Integral 

Absolute  Error 

Relative  Error 

5 

0.246599205561 149E-01 

0.13802E-01 

0.35884E+00 

■a 

0.384268064253340E-01 

0.34732E-04 

0.90303E-03 

15 

0.384615378508424E-01 

0.61070E-09 

0.15878E-07 

20 

0.384615384615384E-01 

0.11102E-15 

0.28866E-14 
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Table  18:  Integration  by  the  Gaussian  Quadrature  for  the  Products  of 
Polynomials  and  Fractional  Power  s(x)  =  x* 
in  TABLE  7 

/(x)dx  =  ^  u></(xi) 


/(x)  =  sin(lOx) 


fix)dx  =  0.183907152907645E+00 

N 

Computed  Integral 

Absolute  Error 

Relative  Error 

5 

0.2023540383881 79E+00 

0.18447E-01 

0.10031E+00 

10 

0.183906300136278E+00 

0.85277E-06 

0.46370E-05 

15 

0.183907152908747E+00 

0.11020E-11 

0.59920E-11 

20 

0.183907152907645E+00 

0.19429E-15 

0.10565E-14 

fix)  =  x25 

/d  f{x)dx  =  0.384615384615385E-01 


N 

Computed  Integral 

Absolute  Error 

Relative  Error 

5 

0.242788857020719E-01 

0.14183E-01 

0.36875E+00 

10 

0.384243967690490E-01 

0.37142E-04 

0.96568E-03 

15 

0.384615377784612E-01 

0.68308E-09 

0.17760E-07 

20 

0.384615384615385E-01 

0.41633E-16 

0.10825E-14 

/(x)  =  x^®s(x) 

f{x)dx  =  0.380228136882129E-01 


N 

Computed  Integral 

Absolute  Error 

Relative  Error 

5 

0.242788857020719E-01 

0.14183E-01 

0.36875E+00 

10 

0.384243967690490E-01 

0.37142E-04 

0.96568E-03 

15 

0.384615377784612E-01 

0.68308E-09 

rnsmmm 

20 

0.384615384615385E-01 

0.41633E-16 

0.10825E-14 

Table  19:  Integration  by  the  Gaussian  Quadrature  for  the  Products  of 
Polynomials  and  Fractional  Power  s(x)  = 
in  TABLE  8 

/  f{x)dx  =  Yiwif{xi) 

•'0 


/(x)  =  sin(lOx) 


/o'  fix)dx  =  0.183907152907645E+00 

N 

Ck>mputed  Integral 

Absolute  Error 

Relative  Error 

5 

0.185981461566628E+00 

0.20743E-02 

0.11279E-01 

10 

0J83906491673941E+00 

0.66123E-06 

0.35955E-05 

15 

0.183907152908381E+00 

0.73541E-12 

0.39988E-11 

20 

0.183907152907644E+00 

0.80491E-15 

0.43767E-14 

fix)  =  x»5 


/o‘  /(x)dx  =  0.384615384615385E-01 

N 

Computed  Integral 

Absolute  Error 

Relative  Error 

5 

0.218060230285103E-01 

0.16656E-01 

0.43304E+00 

m 

0.384060797524605E-01 

0.55459E-04 

0.14419E-02 

15 

0.384615371307963E-01 

0.13307E-08 

0.34599E-07 

20 

0.38461538461 5383E-01 

0.19429E-15 

0.50515E-14 

fix)  =  x“s(x) 


/o'  /(x)dx  =  0.380228136882129E-01 

N 

Computed  Integral 

Absolute  Error 

Relative  Error 

5 

0.218060230285103E-01 

0.16656E-01 

0.43304E+00 

10 

0.384060797524605E-01 

0.55459E-04 

0.14419E-02 

15 

0.384615371307963E-01 

0.13307E-08 

20 

0.384615384615383E'01 

0.19429E-15 

0.50515E-14 
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Table  20:  Integration  by  the  Gaussian  Quadratiire  for  the  Products  of 
Polynoniials  and  Fractional  Power  s{x)  = 
in  TABLE  9 


I  /(i)dx  =  53tui/(i,) 

»=i 


/(x)  =  sin(lOx) 


fix)dx  =  0.182 

1907152907645E+00 

N 

Computed  Integral 

Absolute  Error 

Relative  Error 

5 

0.181420301850517E+00 

0.24869E-02 

0.13522E-01 

10 

0.183906605077266E+00 

0.54783E-06 

0.29788E-05 

15 

0.183907152908166E+00 

0.52114E-12 

0.28337E-11 

20 

0.183907152907645E+00 

0.30531E-15 

0.16601E-14 

fix)  =  x« 


fix)dx  =  0.384615384615385E-01 

N 

Computed  Integral 

Absolute  Error 

Relative  Error 

5 

0.213608283352790E-01 

0.17101E-01 

0.44462E+00 

10 

0.384022552100793E-01 

0.59283E-04 

0.15414E-02 

15 

0.38461536975381 7E-01 

0.14862E-08 

0.38640E-07 

20 

0.38461538461 5384E-01 

0.69389E-16 

0.18041E-14 

fix)  =  x*®s(x) 


fi  fix)dx  =  0.389105058365759E-01 

N 

Computed  Integral 

Absolute  EIrror 

Relative  Error 

5 

0.213608283352790E-01 

0.17101E-01 

0.44462E+00 

10 

0.384022552100793E-01 

0.59283E-04 

0.15414E-02 

15 

0.38461536975381 7E-01 

0.14862E-08 

0.38640E-07 

20 

0.38461538461 5384E-01 

0.69389E-16 

0.18041E-14 
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Table  21:  Integration  by  the  Gaussian  Quadrature  for  the  Products  of 
Polynomials  and  Fractional  Power  s(x)  =  x~^ 
in  TABLE  10 

[  f{x)dx  =  '^Wif{xi) 


f{x)  =  sin(lOx) 


/(x)dx  =  0.183907152907645E-fOO 

N 

Computed  Integral 

Absolute  Error 

Relative  Error 

5 

0.170313857697349E+00 

0.13593E-01 

0.73914E-01 

10 

0.183906941803763E+00 

0.21110E-06 

0.11479E-05 

15 

0.183907152907514E+00 

0.13078E-12 

0.71114E-12 

20 

0.183907152907644E+00 

0.14155E-14 

0.76970E-14 

fix)  = 


/o*  /(x)dx  =  0.384615384615385E-01 

N 

Computed  Integral 

Absolute  Error 

Relative  Error 

5 

0.204389144887852E-01 

0.18023E-01 

0.46859E+00 

10 

0.383937885473486E-01 

0.67750E-04 

0.17615E-02 

15 

0.384615366085782E-01 

0.18530E-08 

0.481 77E-07 

20 

0.384615384615386E-01 

0.11102E-15 

0.28866E-14 

fix)  =  x“s(x) 


/o*  fix)dx  =  0.392156862745098E-01 

N 

Computed  Integral 

Absolute  Error 

Relative  Error 

5 

0.214450816%4760E-01 

0.17771E-01 

0.45315E+00 

10 

0.391572846318390E-01 

0.58402E-04 

0.14892E-02 

15 

0.392156850446440E-01 

0.12299E-08 

0.31362E-07 

20 

0.392156862745100E-01 

0.16653E-15 

0.42466E-14 
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Table  22:  Integration  by  the  Gaussian  Quadrature  for  the  Products  of 
Polynomials  and  Fractional  Power  s(x)  =  x~^ 
in  TABLE  11 


[  f{x)dx  =  '^Wif{xi) 
•*! 


/(x)  =  sin(lOx) 


/o‘  /(x)dx  =  0.183907152907645E+00 

N 

Computed  Integral 

Absolute  Error 

Relative  Error 

5 

0.156179751994470E+00 

0.27727E-01 

0.15077E+00 

10 

0.183907466630775E+00 

0.31372E-06 

0.17059E-05 

15 

0.183907152906471E+00 

0.11742E-11 

0.63847E-11 

20 

0.183907152907646E+00 

0.10825E-14 

0.58859E-14 

m = 


f{x)dx  =  0.384615384615385E-01 

N 

Computed  Integral 

Absolute  Error 

Relative  Error 

5 

0.194693376078146E-01 

0.18992E-01 

0.49380E+00 

10 

0.3838406968 !1997E-01 

0.77469E-04 

0.20142E-02 

15 

0.384615361512299E-01 

0.23103E-08 

0.60068E-07 

20 

0.384615384615379E-01 

0.52042E-15 

0.13531E-13 

/(x)  =  x“s(x) 


si  f{x)dx  =  0.395256916996047E-01 

N 

Computed  Integral 

Absolute  Error 

Relative  Error 

5 

0.194693376078146E-01 

0.49380E+00 

10 

0.383840696841997E-01 

0.77469E-04 

0.20142E-02 

15 

0.384615361512299E-01 

0.23103E-08 

0.60068E-07 

20 

0.384615384615379E-01 

0.52042E-15 

0.13531E-13 

p 
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